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Abstract 

This tutorial review gives an overview of the transverse Anderson localiza¬ 
tion of light in one and two transverse dimensions. A pedagogical approach is 
followed throughout the presentation, where many aspects of localization are il¬ 
lustrated by means of a few simple models. The tutorial starts with some basic 
aspects of random matrix theory, and light propagation through and reflection 
from a random stack of dielectric slabs. Transverse Anderson localization of light 
in one- and two-dimensional coupled waveguide arrays is subsequently estab¬ 
lished and discussed. Recent experimental observations of localization and image 
transport in disordered optical fibers are discussed. More advanced topics, such 
as hyper-transport in longitudinally varying disordered waveguides, the impact 
of nonlinearity, and propagation of partially coherent and quantum light, are also 
examined. 
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1. Introduction 


Anderson localization is the absence of diffusive wave transport in highly disordered 
scattering media [1-4]. It was first introduced in a pioneering theoretical study in 1958 
by Philip Warren Anderson [1], who investigated the behavior of spin diffusion and 
electronic conduction in random lattices. It took more than ten years for the scientific 
community to recognize the importance of Anderson’s work. However, it has remained 
at the forefront of physics research since 1968. There are still many uncertainties and 
unanswered questions in the linear and nonlinear behavior of disordered systems in 
various dimensions. 

The model that Anderson studied involved an electron on a potential lattice with a 
random spread in the energies of the sites caused by a source of disorder. The electron 
was allowed to hop between sites via nearest neighbor potential coupling terms. An¬ 
derson showed that the wavefunction of electron localizes to only few sites at all times, 
provided that the amount of randomness is sufficiently large. 

It did not take long for Anderson and others to realize that the novel localization phe¬ 
nomenon was due to the wave nature of the quantum mechanical electrons scattering 
in a disordered potential, and that similar behavior should also be observed in other co¬ 
herent wave systems, including classical ones [4-8]. The disorder-induced localization 
in electronic systems was shown to be inhibited by thermal fluctuations and nonlin¬ 
ear effects; therefore, it was encouraging to find other avenues in which the disorder- 
induced Anderson localization could be observed. Subsequently, localization was stud¬ 
ied in various classical wave systems including acoustics, elastics, electromagnetics, 
optics [4-11], and various quantum optical systems, such as atomic lattices [12] and 
propagating photons [13-16]. 

Optical systems have played a unique role in the fundamental understanding and ex¬ 
perimental observation of Anderson localization. Optical phenomena are easy to “vi¬ 
sualize,” and there are many advanced tools and techniques in optics that can be used 
to study the physics of localization. Optical studies of Anderson localization can often 
be done with tools that are widely accessible and can be performed in a single labo¬ 
ratory. In addition, Anderson localization already has device-level applications in the 
optics [17-19], and optics can “illuminate” the path to localization-based devices in 
other disordered classical and quantum wave systems. 

It has been shown that coherent waves in one-dimensional (ID) and and two- 
dimensional (2D) unbounded disordered systems are always localized [20]. For 
bounded ID and 2D systems, if the sample size is considerably larger than the localiza¬ 
tion radius, the boundary effects are minimal and can often be ignored [21,22]. How¬ 
ever, in three-dimensional (3D) coherent wave systems, the scattering strength needs 
to be larger than a threshold value for the localization to happen [23]. The scattering 
strength is characterized by the wave scattering transport length /* (shorter l* means 
stronger scattering), and the Ioffe-Regel condition [24] states that in order to observe 
Anderson localization, the disorder must be strong enough that the wave scattering 
transport length becomes on the order of the wavelength. The Ioffe-Regel condition is 
often cast in the form of kl* ~ I, where k is the effective wavevector in the medium. 
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It is notoriously difficult to satisfy in 3D disordered-media. For example, for the 
optical field to localize in 3D, very large refractive index contrasts are required that 
are not generally available in low-loss optical materials [7]. “The fact that Anderson 
localization is hard to achieve in 3D optical systems may be a blessing in disguise; 
otherwise, no sunlight would reach the earth on highly cloudy days” [25]. In order 
to observe Anderson localization of light, strongly scattering materials at optical and 
near infrared frequencies such as TiCF, GaAs, GaP, Si, and Ge nanoparticles can be 
used. Careful measurements are required because bulk absorption can easily lead to 
experimental signatures similar to Anderson localization [26-28]. 

Unlike 3D lightwave systems, in which observation of localization is prohibitively 
difficult, observation of Anderson localization in quasi-2D and -ID optical systems 
(transverse Anderson localization) is readily possible, as was first shown by Abdullaev 
et al. [29] and De Raedt el al. [30]. There have since been many reports on the observa¬ 
tion of transverse Anderson localization of light in ID and 2D, which is the main focus 
of this tutorial review, as well. Transverse Anderson localization is attractive because 
of its relative simplicity, ease of experimentation, and the rewarding physical insights it 
brings about on many fronts. Moreover, the longitudinal coordinate along the direction 
of propagation plays the role of time in a 2D disordered system; therefore, controlled 
temporal variations can also be studied in these systems. 

There are many excellent reviews that cover various aspects of wave propagation 
in disordered systems as well as Anderson localization [3,7,26,31-37], This tutorial 
review is neither intended to be comprehensive, nor is it intended to explore each area 
it covers in great depth. Rather, the intent is to provide a pedestrian and intuitive ap¬ 
proach to Anderson localization, mainly focused on the transverse localization of light. 
The coverage of topics is inevitably slanted toward those of particular interest to the au¬ 
thor. The pedagogical approach is intended to benefit both newcomers to this rewarding 
research field, as well as outsiders who are interested to learn about Anderson local¬ 
ization. The author seeks forgiveness from those whose work is not mentioned here, as 
well as for any technical errors or omissions. 

2. A random matrix example 

A good way to build an intuition about the relationship between randomness and lo¬ 
calization is to use random matrices [36,38]. The following example shows that the 
extended eigenvectors of an ordered matrix become very localized when some ran¬ 
domness is added to the elements of a matrix. 

Consider a symmetric tridiagonal N x N matrix M defined as 

M i.i = 1, M,y + i = = 0.1, (1) 

for all possible values of i. For definiteness in this numerical example, we consider 
N = 200. Matrix M has N real eigenvectors, and each eigenvector is an /V-clemcnt 
vector. We identify the ith eigenvector as V (j K and represents its /th element. 

In Figure 1, a few eigenvectors V ' 1 are plotted as a function of their element number 
j. In Figure 1(a), Vy' 1 , , V'p , and are plotted and are all oscillatory functions 
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of j, where the zth eigenvector is identified by / + I oscillation nodes. V^"°° ; is also 
plotted in Figure 1(b). The relevant observation here is that all of these eigenvectors 
are extended over the entire element-position domain. 

We now would like to show that randomness can localize the eigenvectors of M over 
the element-position domain. For this part, lets preserve the tridiagonal character of M 
as presented in Eq. 1, but add a small random number to each off-diagonal element. 
The new, slightly randomized M is formally defined as 

M,y = 1, M,y + 1 =M ; - + i ji - = 0.1 + r i , Ti £ unif[—0.01,0.01], (2) 

where r,- is a number randomly selected from a real uniform distribution in the range 
[—0.01,0.01]. Similar to the case of the ordered M in Figure 1, and 

V^ are plotted in Figure 2(a) and ¥^° n is plotted in Figure 2(b). This time, all of 
these eigenvectors appear to be localized over the element-position domain. 

Not only does the randomness result in localized eigenvectors, but also the relative 
strength of randomness compared with the average off-diagonal values determines the 
size of the localization (localization length). This can be verified by increasing the 
range of the random numbers r,- to [—0.05,0.05]. The new, strongly randomized M is 
formally defined as 

M,y = l, = Mj+q,- = 0.1 + r,-, r, £ unif[—0.05,0.05]. (3) 

Again, the same eigenvectors are plotted in Figures 3(a) and 3(b). The eigenvectors are 
strongly localized over the element-position domain because of the strong randomness 
in the off-diagonal elements of M. 

The impact of disorder on the distribution of the width of the eigenvectors on the 
element-position domain can be easily visualized in a histogram. In Figure 4, the dis¬ 
tribution of the width of the eigenvectors is plotted for the weak disorder of Eq. 2 and 
the strong disorder of Eq. 3. The width is calculated using the second moment method. 




FIGURE 1. This figure shows that the eigenvectors of the ordered matrix M de¬ 
fined in Eq. 1 are extended over the entire element-position domain. Eigenvectors 
andVy 4 ' are plotted in (a); andV^ 200 * is plotted in (b) as a function 
of the element position j. 
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FIGURE 2. Similar to Figure 1, except the matrix M is slightly randomized ac¬ 
cording to Eq. 2 and the eigenvectors are localized. 


which is the standard deviation around the mean position calculated using the absolute- 
value-squared of the eigenvectors. The width a,- for the eigenvector VW is given by 



Each probability distribution is presented in a histogram and is the result of averaging 
over 100 independent random simulations. From the distributions in Figure 4, it is clear 
that stronger disorder results in stronger localization of the eigenvectors; moreover, it 
is clear that such statements can only be made in a statistical sense. Therefore, in the 
case of strong disorder, although the majority of the eigenvectors are more localized, 
a minority of the eigenvectors may actually be less localized than those of the weak 
disorder. 

The exercise presented here shows how off-diagonal disorder results in localization. 
Similar localization behavior can be observed for diagonal disorder, where only the 
diagonal elements of M are randomized, and also for mixed diagonal and off-diagonal 
disorder. 
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FIGURE 3. Similar to Figure 1 and Figure 2, except the matrix M is strongly 
randomized according to Eq. 3 and the eigenvectors are strongly localized. 
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FIGURE 4. The probability distribution of the eigenvector widths for two separate 
cases of the weak disorder related to Eq. 2, and the strong disorder related to Eq. 3. 



FIGURE 5. (a) A periodic array of two different dielectric materials identified 
with refractive indexes n\ and rci is shown, (b) is similar, except the thickness of 
the layers is chosen randomly. 


Highlights: 

• The eigenvectors of random matrices can be localized in the element 
position space. 

• Some modes are very narrow and some are wide, and localization is 
only meaningful in a statistical sense. 

• By calculating the width of the eigenvectors of a large ensemble of 
random matrices, it is possible to calculate the probability distribution 
for the width of the eigenvectors. 

• A stronger level of randomness shifts the eigenvector-width probabil¬ 
ity distribution to smaller width values, hence a stronger localization. 


3. Normal transmission through a random stack of dielectrics 

Another interesting example that links randomness to localization is the problem of 
light transmission through a random stack of dielectrics. The normalized transmission 
through a stack of dielectrics is shown in Figure 5. The light, which is incident from 
the left, is partially reflected from the stack, while the rest is transmitted through the 
stack. The dielectrics are assumed to be lossless. Figure 5(a) shows a periodic array 
of two different dielectric materials identified with refractive indexes n\ and n 2 , more 
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commonly referred to as a Bragg grating. Figure 5(b) is similar, except the thickness of 
the layers is chosen randomly. 

In Figure 6(a), the relative optical power transmission is plotted as a function of 
the normalized frequency. The thickness of each layer is A (identical for all layers), 
and ko = 2n /A is the wavevector, where A is the optical wavelength in vacuum. The 
stack is made of 200 air-glass layers (400 layers total), where n\ = 1.5 and n 2 = 1.0 is 
assumed. The usual bandgaps and bandpasses are observed in the transmission plot of 
this periodic Bragg grating. 

The setup for Figure 6(b) is identical to that of Figure 6(a), except the thickness of 
each layer is chosen from a uniform random distribution in the range [0,2A]. Unlike 
the periodic Bragg grating, transmission through the random stack is nearly zero except 
for very small values of the normalized frequency kc,A. Therefore, the random stack 
behaves like a nearly perfect mirror beyond a certain frequency. 




FIGURE 6. The optical transmission through a stack of 200 glass-air layers is 
plotted versus k ()A for (a) a periodic layer thickness and (b) random layer thick¬ 
ness. A is the periodicity in case (a), while the thickness of each layer in case (b) 
is chosen from a uniform random distribution in the range [0,2A]. 


The calculation of the optical transmission in Figure 6 is carried out using the trans¬ 
mission matrix M defined in Eq. 5 related to Figure 7 : 



(n 2 +n l )e' <p 
(n 2 — n\)e~ 1<? 


(n 2 — n i)e l<p 
(ri 2 +ni)e~ l<p 


(p = n\k()d. 


(5) 


In Figure 6, the M matrix can be used to relate the right- and left-moving components 
of the optical field in dielectric n 2 to those in dielectric n\ according to 


u? 


= M. 


ur 


( 6 ) 


The total transmission and reflection can be calculated by cascading all the M matrices 
for different random layers of the dielectric stack, as discussed in detail in Ref. [39]. 

It is instructive to repeat the exercise of Figure 6 with a larger number of layers. 
In Figure 8, the number of layers is increased to 10,000 in each case, and the relative 






























FIGURE 7. A sketch of the field amplitudes, refractive indexes, and geometry 
related to Eq. 5 and Eq. 6. 


optical power transmission is plotted in dB units. The logarithmic scaling on the verti¬ 
cal axis shows that the bandgap frequency ranges in the ordered stack of Figure 8(a), 
although not perfect due to the finite number of layers, attenuate the transmitted light 
by more than 20,000 dB, making the stack a nearly perfect mirror over this frequency 
range. The random stack in Figure 8(b), except for low frequency values, is also a very 
good reflector. Although the attenuation in trans mi ssion at 3,000 dB or more is not as 
good as the bandgap region in the periodic case, it is very broadband and is not limited 
to bandgap range. 

The observations in Figure 8 bring about a somewhat philosophical issue. In prac¬ 
tice, any stack of optical dielectrics has some inevitable randomness; therefore, given 
a sufficiently large number of layers, it acts as a perfect mirror beyond a certain fre¬ 
quency. Of course, at very low frequencies, the wavelength is so large that it hops over 
the stack and results in large transmission. Therefore, in practice, any stack of opti¬ 
cal dielectrics with a sufficiently large number of layers is a practically a broadband 
bandgap structure. The periodic structure analyzed in Figure 8(a) can be viewed as the 
limiting case of the random structure, when the randomness goes to zero. When this 
limit is taken, the bandpasses emerge out of the broadband bandgap of the more gen¬ 
eral disordered stack. Therefore, the philosophical point of view is that the magic of 
a periodic structure is not in its bandgap, because the bandgap comes naturally in any 
practical ID stack; rather, it is the emergence of the bandpasses that makes periodic 


10000 glass-air layers, ordered 



10000 glass-air layers, random 



FIGURE 8. Same as Figure 6, except with 10,000 layers. 
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structure so special. While this may seem like an inconsequential philosophical issue, 
it is necessary to alert students who learn about coherent waves in periodic systems in 
optics or condensed matter physics to ensure that they develop proper early intuition 
about the true implications of periodicity. The bandpasses are at least as glorious as 
the bandgaps. 

The low-frequency behavior, as well as the periodic “very weak” resurrection of 
transmission in Figure 8(b), can be intuitively understood based on the work of Berry 
et al. [40]. They have shown that the optical transmission through a random stack of 
N bi-layered transparent plates, each contributing a random phase 0 belonging to a 
uniform distribution 0 e unif[0,27r], is given by 

T M = exp( — 2NIog(l/0), (7) 

Using N = 10,000, n\ = 1.5, and «2 = 1.0, one obtains X 2 N ~ —3546 dB, which is in 
agreement with the numerical simulation in Figure 8(b). In the simulations presented 
in Figures 6 and 8, it was assumed that {d\ = do) G unif[0,2A]; therefore, 

0 < 0 = (nikod + mkod) < 2(n\ + n 2 )k()A. (8) 

In order to obtain the near uniform distribution of 0 e unif[0, 2n , the upper bound in 
Eq. 8 can be set to 2n, where one obtains 

k 0 = k 0 A^n/(ni+n 2 ), (9) 

as the condition for the disorder-induced localization, exerting its full power and agree¬ 
ing with the assumptions of Ref. [40]. It should be noted that this value is very close 
to the value of koA in Figure 8(b), beyond which the attenuation is strong. Therefore, 
the main reason that the attenuation is not strong at low /cqA is that the random phase 
0 does not cover the entire range of [0,2 jt], For koA > ko, although the range of the 
phase 0 covers the entire [0,2 jt], it partially folds over and makes the distribution of 0 
non-uniform over [0,2 jt]. 

For example, consider the case of 0 £ unif[0.3 tt] . 0 is a phase variable for which 
[271,371] = [0,7r]; therefore, 0 covers the entire [0, 2 tt] range but with twice the proba¬ 
bility in the [0,7r] range compared with \n. 2n range and is no longer uniformly dis¬ 
tributed over [0, 2 tt] . This is the main reason behind the partial resurrections of trans¬ 
mission observed in Figure 8(b) for k$A > ko- 

A final and important point is that the exponential decay of the optical amplitude in 
the disordered stack is a manifestation of the “coherent” superposition of partial waves 
that reflect and transmit at each boundary. The fact that the naive ray theory, which is 
based on incoherent transmissions and reflections, fails to predict the Anderson local¬ 
ization behavior observed in these examples= defies the intuition at first glance. The 
interested reader can consult Ref. [40] for a more detailed discussion on these points. 
A simple and elegant experiment using a laser pointer and a stack of overhead trans¬ 
parencies is also presented in Ref. [40], which can be excellent demonstration in an 
undergraduate laboratory in optics. 
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Highlights: 


• Optical transmission through a stack of dielectric layers with differ¬ 
ing refractive indexes and random thicknesses drops, on the average, 
exponentially with the number of layers. 

• This exponential drop is observed for all wavelengths that are compa¬ 
rable or shorter than the mean thickness of the dielectric layers. The 
random stack behaves like an ultra-broadband mirror. 

• The optimum localization is obtained when the cumulative optical 
phase in each layer is from a uniformly distributed random number 
in the range [0,2 tt]. If the domain is smaller, or even larger such as 
[0,3 k\, the localization is weaker. 


4. Oblique transmission through a random stack of dielectrics 

The example presented in section 3 can be readily generalized to the case of light 
incident at an angle on the dielectric stack, as illustrated in Figure 9. The calculation 
can be carried out using a simple generalization of the transmission matrix M defined 
in Eq. 5 to the case of incidence at an angle (see for example Ref. [39]). 

In Figure 10, the relative optical power transmission is plotted as a function of the 
normalized frequency for oblique incidence at a 45° angle (as measured in the air layer) 
for the TE (transverse electric) polarization of light. Figure 10 should be compared with 
Figure 8, which was for the case of normal incidence. 

In the case of the periodic dielectric stack in Figure 10(a), the bandgaps are shifted 
compared with normal incidence in Figure 8(a) and are wider and deeper. For the ran¬ 
dom dielectric stack in Figure 10(b), the transmission is lower over the entire frequency 
band compared with Figure 8(b); however, besides the lower value and frequency shifts 
of the features in transmission, no substantial qualitative difference exists between 
oblique incidence at a 45° angle and normal incidence. 

The case of a TM (transverse magnetic) polarization is explored in Figure 11, where 
the periodic-stack bandgaps in Figure 11(a) are substantially narrower and more shal¬ 
low compared with the case of a TE polarization at 45° angle in Figure 10(a) and 
normal incidence in Figure 8(a). Similarly, the transmission of the TM polarization at 
45° angle in Figure 11(b) through a random stack is stronger in the case of a TM po- 



FIGURE 9. (a) Transmission and reflection of light incident on a periodic stack 
of dielectrics, (b) Same, except the thickness of each dielectric layer is randomly 
selected. 
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FIGURE 10. Relative optical power transmission is plotted as a function of the 
normalized frequency for oblique incidence at a 45° angle for the TE polarization 
for (a) a periodic dielectric stack and (b) a random dielectric stack. 


larization compared with the TE polarization in Figure 10(b) and the case of normal 
incidence in Figure 8(b). 

For the larger incidence angle of 85°, the difference between transmission through a 
periodic dielectric stack and a random dielectric stack is more pronounced; and so is the 
difference between TE and TM poalizations. The bandgaps are wider and deeper (note 
the 10 4 dB label in the vertical scale) for the TE polarization incident on a periodic 
dielectric stack in Figure 12(a) compared with the case of 45° and the case of normal 
incidence. The gaps are also shifted in normalized frequency and are more widely sep¬ 
arated. Similar behavior is observed for the case of the TE polarization incident on a 
random dielectric stack in Figure 12(b), where the full strength of the localization be¬ 
gins at a higher normalized frequency than what is shown in Figure 12(b); this behavior 
agrees with the general shift of the spectral features to the higher normalized frequency 
observed in the case of a periodic dielectric stack in Figure 12(a). 

For the TM polarization incident on a periodic dielectric stack in Figure 1 3(a), many 
namow and closely separated gaps appear in the studied normalized frequency range. 
The gaps are deeper compared with the case of 45° and the case of normal incidence. 
Similarly, the spectral features shift to lower values in the case of the TM polarization 


10000 glass-air layers, ordered 



10000 glass-air layers, random 



0 1 2 3 4 5 

k 0 A 


FIGURE 11. Relative optical power transmission is plotted as a function of the 
normalized frequency for oblique incidence at a 45° angle for the TM polarization 
for (a) a periodic dielectric stack and (b) a random dielectric stack. 


12 





















10000 glass-air layers, ordered 



10000 glass-air layers, random 



FIGURE 12. Relative optical power transmission is plotted as a function of the 
normalized frequency for oblique incidence at an 85° angle for the TE polarization 
for a (a) periodic dielectric stack and a (b) random dielectric stack. 


incident on a random dielectric stack in Figure 13(b); therefore, the full strength of the 
localization begins at a lower normalized frequency compared with TE polarization 
and smaller incidence angles (for both polarizations). 

TM polarization and Brewster’s angle 

In the above examples, it is observed that the TM polarization behaves differently from 
the TE polarization. The root of this distinction is in the different boundary conditions 
for the electric field in the TE and TM polarizations at the interface of each layer, and 
the difference is most apparent at Brewster’s angle, while the TM polarization shows 
perfect transmission. In Figure 14, the relative optical power transmission is plotted as 
a function of the incidence angle of the TM polarized light on a random dielectric stack 
for (a) IcqA = 1, and (b) kg A = 5. Perfect transmission is observed at Brewster’s angle 
6b = 56.3° through 10,000 glass-air layers of random thickness. Of course, like other 
simulations so far, intrinsic attenuation of glass is neglected in these simulations. The 
transmission is lower (localization is stronger) for k^A = 5 in Figure 14(b) compared 
with the case of k^A = 1 in Figure 14(a), as expected from previous simulations and 
related arguments. 

10000 glass-air layers, ordered 10000 glass-air layers, random 




FIGURE 13. Relative optical power transmission is plotted as a function of the 
normalized frequency for oblique incidence at an 85° angle for the TM polariza¬ 
tion for (a) a periodic dielectric stack and (b) a random dielectric stack. 
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10000 glass-air layers, k 0 A =1 10000 glass-air layers, k 0 A =5 




FIGURE 14. Relative optical power transmission is plotted as a function of the 
incidence angle of TM polarized light on a random dielectric stack for (a) koA = 1, 
and (b) koA = 5. Perfect transmission is observed at Brewster’s angle 0b = 56.3°. 


For comparison, in Figure 15, the relative optical power transmission is plotted as a 
function of the incidence angle of TE polarized light on a random dielectric stack for 
(a) ko A = 1, and (b) k$A = 5. Compared with Figure 14, the absence of perfect trans¬ 
mission at Brewster’s angle is notable. It can also be concluded that below Brewster’s 
angle, transmission is weaker in the TM polarization compared with the TE. This is 
expected, because at normal incidence TE and TM polarizations have the same trans¬ 
mission; increasing the angle results in a monotonic increase in transmission for the 
TM polarization all the way up to 100% at Brewster’s angle, while the trans mi ssion for 
the TE polarization monotonically decreases. Beyond Brewster’s angle, the TE polar¬ 
ization continues its slow decline; however, the TM polarization goes through a steep 
decline. Depending on the value of k$A, at large angles, the TM polarization will have 
a lower transmission compared with TE (compare Figure 14(a) with Figure 15(a)), 
or vice versa (compare Figure 14(b) with Figure 15(b)). The interested reader is en¬ 
couraged to consult a detailed account of the polarization dependence of an obliquely 
incident light on layered media in Ref. [41], 
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FIGURE 15. Relative optical power transmission is plotted as a function of the 
incidence angle of TE polarized light on a random dielectric stack for (a) kg A = 1, 
and (b) kg A = 5. 
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Random dielectric stack as a broadband waveguide 

Based on the detailed discussions of the oblique incidence e.g., the results in Figure 14 
and Figure 15, it is conceivable to make a ID waveguide in the form shown in Fig¬ 
ure 16. For a periodic dielectric stack, a waveguide can efficiently operate over a nar¬ 
row range of wavelengths in the bandgaps corresponding to the specific incident wave¬ 
length [42^44]. The operating bandwidth can be increased by considering a chirped 
grating as shown in Ref. [45]. Alternatively, a design based on a random dielectric 
stack is very broadband; the trade off is that the leakage is slightly higher than that of a 
bandgap design for a periodic dielectric stack. Therefore, if ultra-broadband operation 
is desired, a random stack is likely the solution. 



FIGURE 16. A ID dielectric waveguide using (a) a periodic dielectric stack, 
which can have a very low leakage over a narrow range of wavelengths in the 
bandgaps corresponding to the specific incident wavelength and incident angle; 
and (b) a random dielectric stack, which can operate over a broad range of wave¬ 
lengths and angles, but with slightly more leakage compared with a periodic 
waveguide of the same number of layers optimized to operate in the center of 
the bandgap. 
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Highlights: 


• Considerable differences are observed between the transmission of 
the TE and TM polarizations, when a random dielectric stack is il¬ 
luminated at an angle. This hints at the possibility that the vectorial 
nature of the electromagnetic field may play an important role in set¬ 
ting the localization behavior in certain situations. 

• The strongest difference between the TE and TM polarizations can 
be observed for oblique incidence at Brewster’s angle, at which the 
TM polarization is entirely transmitted through the random dielectric 
stack. 

• Below Brewster’s angle, the TE polarization usually has a lower 
transmission; however, above Brewster’s angle, the situation depends 
on the ratio of the size of the optical wavelength to the average lattice 
size. 

• The concept of a broadband waveguide is introduced, where the ran¬ 
dom dielectric stacks are used as broadband reflecting walls for the 
waveguide. 


5. Transverse Anderson localization of light in one transverse dimension 

Consider the one dimensional array of N identical single-mode optical fibers sketched 
in Figure 17. The propagation constant of each fiber is /3o and the direction of the 
propagation of the optical wave is assumed to be into the page. Each fiber is weakly 
coupled to its nearest neighbor, and the strength of the coupling is determined by the 
separation between the fibers. The coupled mode equations for the propagation of the 
optical field through this optical fiber array can be expressed as [39] 



( 10 ) 


A j is the amplitude of the optical field in the jth fiber, c t (cj ) is the coupling strength of 
the /tli fiber to its right neighbor A j + \ (left neighbor Aj _ i), and we assume symmetric 
coupling, so cj = cj +l . Of course, we need to assume cj - = 0 and = 0, because the 
1st (Nth) fiber does not have a neighbor to its left (right). 


j=l j= 2 j= 3 j= 4 j= 5 


j=N-l 




FIGURE 17. A one dimensional array of N identical single-mode optical fibers, 
referred to as a coupled waveguide array. 

Here, we consider the simple case, where the input light is only coupled to the middle 
fiber at z = 0. As the optical field propagates through the middle fiber, it couples to its 
neighboring fibers; those will couple to their neighbors as well, and this cascading 
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event results in a discrete diffraction pattern. This scenario is shown in Figure 18(a), 
where the intensity of the propagating light is plotted as a function of the propagation 
distance and the waveguide number. The parameters used in this simulation are: /3o = 6, 
Co = 0.01, N = 201, and 0 < z < 5000. The input boundary condition is set to Aj{z = 
0) = 5/joi, as mentioned above. 

The situation is quite different if the coupling between the waveguides is random¬ 
ized. For example, consider perturbing the array slightly, so that the average coupling 
remains co, but individual couplings vary with a uniform distribution according to 
c+ = co + rj, where rj E unif[—0.006,0.006]. In Figure 18(b), we observe that the light 
in the middle waveguide cannot spread as efficiently to its neighbors and the diffraction 
slows down as the beam propagates farther in z. This behavior is certainly induced by 
the disorder and randomness introduced in the cross coupling of fibers. Therefore, it is 
reasonable to expect that more randomness should result in slower diffraction. This is 
verified in Figure 18(c) by repeating the experiment for rj E uniff -0.01.0.01. 

Let's take a closer look at the optical field pattern in Figure 18(c). The propagating 
light initially follows a diffracting pattern, but after a certain propagation distance the 
diffraction halts and the total width of the beam remains more or less the same for the 
rest of the propagation. It appears that the fibers near the edges of the array always 
remain dark. This behavior is referred to as transverse Anderson localization. In fact, 
similar behavior can be observed in Figure 18(b) after a sufficiently long propagation 
distance; however, the transition to a stable width occurs at a longer propagation dis¬ 
tance and the beam width is larger because of the smaller amount of disorder. 

The light propagation in a one-dimensional disordered coupled fiber array is trans¬ 
versely localized for any amount of disorder; however, if the disorder is too small, the 
stable localized beam width may be larger than the transverse size of the structure and 
the localization effect cannot be observed in practice. 

The coupled fiber array of Figure 17 described by Eq. 10 can be analyzed alterna¬ 
tively in the language of normal modes. The advantage of the normal mode description 
is that the equations describing light propagation in a disordered coupled fiber array can 



propagation distance (x 100) 


FIGURE 18. Propagation through a waveguide coupled array for the case of (a) 
a periodic array, (b) a disordered array, and (c) a highly disordered array. Higher 
level of disorder results in a more localized propagation. 
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be directly mapped to a random matrix discussed in section 2: the transverse Anderson 
localization is nothing but the localization of the eigenvectors in the element-position 
domain index j of A,-. In order to see this, we define the vector A as 

A = (Ai,A 2 ,A3,-• • ,Ajv-i,Ajv), (11) 

and rewrite Eq. 10 

i-^-A + B.A = 0, (12) 

oz 

where all diagonal elements of the symmetric tridiagonal matrix B are equal to /3o and 
the off-diagonal elements are cj = cq + rj. The boundary condition is set to A j(z = 
0) = 5/.ioi, because the input light is only coupled to the middle fiber. 

The real tridiagonal matrix B is symmetric; therefore, ah of its eigenvectors are real 
and mutually orthogonal and its eigenvalues are real. We identify the ith eigenvalue 
as /3, and the corresponding eigenvector as V-, and V'- ,: represents the /th element 
of the VW eigenvector. Furthermore, without loss of generality, we assume that the 
eigenvectors are properly scaled such that they are orthonormal, i.e., vW.yO') = $‘i. 
Using the eigenvectors, we can construct a real orthogonal matrix Q such that Qjj = 
\ where Q.O 7 ' = I, and cany out an orthogonal transformation to rewrite Eq. 12 as 

/—A + B.A = 0, where A = Q.A, E = Q.B.Q T . (13) 

oz 

The orthogonal transformation of the matrix B results in a diagonal matrix whose ele¬ 
ments are the eigenvalues of B, i.e., Mjj = /3,c),y. 

Using the above information, we can rewrite Eq. 13 as 

i^-Aj(z) + PjAj{z) =0, j = 1,- ■ ,N, (14) 

where we have defined 

A = (Ai,A 2 ,A 3 ,--- ,A n _i,A n ). (15) 

The solutions to the N independent first-order differential equations 14 for the unidi¬ 
rectional propagation of light can be simply written as 

Aj(z)=Aj(0)exp[-iPjz], j = (16) 

Using the orthogonal transformation and a few lines of simple algebra, Eq. 16 can be 
expressed as 


II 

N 

Y, QkJQk,i A i (0)exp[-ij3*z], j = 1 , •• 

■,N. 

(17a) 

k= 1 

i=i 



N 




= L 

QkjQk,m exp[—ijSjtz], 


(17b) 

k= i 




ii 

vfvSi exp [-ipiczj. 


(17c) 
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Equation 17c brings the problem to its final form, using which we can now discuss 
the transverse Anderson localization in terms of the random matrix B in the language 
of Section 2. Recall that we showed in Figure 2 and the corresponding discussions that 
the eigenvectors of a symmetric tridiagonal matrix with random off-diagonal elements 
arc localized over the element-position domain. Therefore, only a few elements of each 
in Eq. 17c are non-zero. We also observed in Figure 2 that the localization of the 
eigenvectors occurs at different locations in the element-position domain; therefore, 
only very few eigenvectors Y- 1 have non-zero elements in the j = 101 position. As 
a result, the sum in Eq. 17c is practically limited to only a few eigenvectors (modes) 
yW that have non-zero j = 101 elements and the other non-zero elements of these few 
eigenvectors are all concentrated around j = 101 due to localization. Therefore, Aj(z) 
will be equal to zero for all z if j is very different from 101. 

From the above discussions we can conclude that for a stronger disorder, the local¬ 
ization becomes stronger as seen in Figure 2, and Aj (z) will remain zero unless j is 
even closer to 101. The picture portrayed here is consistent with Figure 18, especially 
in describing the difference between Figure 18(b) and Figure 18(c). 

Fets recap these observations in the more common modal language. The initial beam 
only couples to those guided modes that arc localized in the vicinity of the center 
fiber. Near the entrance, the optical excitation closely resembles the in-coupling beam, 
as a result of a coherent sum of the excited modes. Each excited mode propagates 
with a different phase velocity determined by the propagation constant (5^, therefore 
the detailed balance between the excitation amplitudes of the guided modes that is 
responsible for the narrow excitation at the entrance Aj(z = 0) = 5/,ioi is broken as the 
relative phases between the modes change when the beam propagates along the fiber 
array. As the beam propagates and the detailed balance is further broken, the beam 
expands; however, it can never expand beyond the size of its constituent modes, which 
are only a few, each being very localized. Therefore, the expansion is eventually halted, 
as observed in Figure 18(b) and Figure 18(c). 

Using the modal language we can explain a few other interesting phenomena. One 
question that is often asked is what happens if the initial excitation is wider than the 
typical width of an individual mode? The answer is that a wider beam excites a few 
more modes; therefore, the propagating light eventually localizes to a slightly wider 
beam. Another interesting fact is that because in Eq. 17c only a few terms in the sum 
play a role, after a sufficiently long propagation distance z r , the beam refocuses back 
to its initial intensity distribution (self-imaging) [46,47], This revival distance can be 
calculated by noting that at z r , (j8* — Pioi )z r = 2 nm^ must be satisfied for all the rel¬ 
evant terms in the sum, where each is an integer. The intensity pattern at z, r will 
be identical to that at z = 0. In practice, if the number of excited modes is sufficiently 
large, the revival distance can be much longer than the waveguide and the refocusing 
is not observed. 

It must be noted that a similar localization behavior is obtained if instead of ran¬ 
domizing the coupling coefficients between the fibers in the array [48], the propagation 
constants of individual fibers are slightly randomized [49]. This was briefly discussed 
in Section 2 as diagonal localization. 
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Mode shapes 

In Figure 19 we show the distribution of the propagation constants of the eigenmodes 
of the coupled waveguide array of Figure 17; i.e., f5/ in Eq. 16, for j = 1, ■ ■ ■ ,201. The 
parameters used for the disorder relate to the case discussed in Figure 18(b), where we 
used /3o = 6, co = 0.01, IV = 201, and ct = CQ + rj, with rj € unif[—0.006,0.006]. The 
calculated modes are numbered from j = 1 to j = N, sorted in a descending order of the 
value of the propagation constant. As expected, the calculated propagation constants 
are in the vicinity of /3o = 6. The range of variation of /j, is determined by the range of rj 
(unif[—0.006,0.006]), while the exact values vary for each different random waveguide. 



FIGURE 19. Propagation constants of a random coupled waveguide array, where 
the modes are numbered from 1 to N = 201, sorted in a descending order of the 
value of the propagation constant. 

In Figure 20 we plot the shape of a few modes of the waveguide. Of course, the 
coupled waveguide array supports a total of N = 201 modes. In Figure 20 each mode is 
normalized to unity for easier comparison. The horizontal axis is the waveguide number 
in Figure 17. Figures 20(a), (b), (c), (d), (e), (f) correspond to mode numbers 1, 25, 
99, 101, 170, 201, respectively. Recall that the mode numbers here correspond to the 
mode numbers used in Figure 19. Therefore, Figure 20(a) is the mode with the largest 
propagation constant in Figure 19, which is the top edge of the propagation constant 
band. The modes that correspond to the region in the vicinity of the top edge of the band 
are highly localized with very few oscillations. As the propagation constant is increased 
e.g., for mode number 25 the mode profile oscillates a few times as in Figure 20(b). The 
modes near the middle of the band in the vicinity of j = 100 generally oscillate rapidly 
and spread over many waveguides, as can be seen for modes number 99 and 101 in 
Figures 20(c) and (d), respectively. These modes become more localized if the disorder 
is increased. Once the mode number gets closer to the bottom edge of the band, the 
oscillations become so rapid that the sign of the mode profile flips between adjacent 
waveguides, and the modes localize again. This behavior can be clearly seen for modes 
number 170 and 201 in Figures 20(e) and (f), respectively. 

The mode width is plotted against the propagation constant of each mode in Fig¬ 
ure 21, averaged over 1000 simulations for (a) rj e unif[—0.006,0.006] and (b) rj e 
unif[—0.01,0.01]. It is clear that the band edges result in the narrowest modes, while 
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waveguide number 


FIGURE 20. Mode profiles of the random coupled waveguide array of Figure 17 
are plotted. Subfigures (a), (b), (c), (d), (e), (f), correspond to mode numbers 1, 25, 
99, 101, 170, 201, respectively. The modes whose propagation constants belong to 
the region near the top edge of the band in Figure 19 are highly localized with no 
or few oscillations as in (a) and (b). The modes with propagation constants near 
the middle of the band spread over many waveguides and oscillate as in (c) and 
(d). The modes with propagation constants near the bottom of the band are highly 
localized and oscillate so rapidly that the sign of the mode profile flips between 
adjacent waveguides as in (c) and (d). 


the modes near the middle of the band are more spread out. Also, the stronger disorder 
in Figure 21 (b) results in smaller mode widths compared with the weaker disorder case 
in Figure 21(a). 




propagation constant propagation constant 

FIGURE 21. The mode width is plotted against the propagation constant 
of each mode. The results are averaged over 1000 simulations for (a) r ; £ 
unif[-0.006,0.006] and (b) rj £ unif[-0.01,0.01]. 


The disorder in the model studied here is off-diagonal. Although the effects of diag¬ 
onal and off-diagonal disorder are generally quite similar, there exist subtle differences 
between the two types of disorder. For example, it has been argued that while all modes 
are exponentially localized for diagonal disorder, this may not be strictly true for off- 
diagonal disorder [50-52]. Specifically, the modes at the center of the band in a ID 
off-diagonally disordered lattice fall off asymptotically as exp(—A \/\j — jo\) instead 
of the simple exponential form exp(— X\j — jo|), where | j — jo represents the distance 
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from the center of the mode at jo. The exp(—Xy/\j — jo\) form is also responsible 
for the divergence of the mode width in the middle of the band in both subfigures of 
Fig. 21 . We emphasize that even in the middle of band, the mode is localized, although 
the weaker localization form results in a divergent mode width. There are also differ¬ 
ences between diagonal and off-diagonal disorder in 2D that will be pointed out later 
in section 6. 

Beam width in position and momentum space 

In Figures 22, 23, 24, 25, we explore in some detail the expansion and localization of 
the optical beam that is coupled to the waveguide array of Figure 17, as a function of 
the longitudinal coordinate. The vertical axis is the beam width, which is determined by 
the second moment method of Eq. 21. The parameters used for the disorder relate to the 
cases already discussed in Figure 18. For these figures we have used j8o = 6, co = 0.01, 
and N = 201. Any figure labeled with (a) corresponds to the disorder-free periodic 
array; labeled with (b) corresponds to the weakly disordered case of cj = co + rj, where 
rj E unif[—0.006,0.006]; and labeled with (c) corresponds to the strongly disordered 
case of rj E uniff—0.01,0.01]. The only minor differences are that 0 < z < 8000, and 
the boundary condition is set to 

Aj(z = 0)=exp[— (j ~^ 2 1)2 ], j = !,■■■ ,N, W 0 = V2. (18) 

We note that the beam width formula of Eq. 21 (presented later in this article) also 
gives the value of y/2 for the beam width at z = 0. 

Figure 22(a) shows the expansion of the beam width in the disorder-free lattice. The 
expansion is ballistic, i.e., the width grows linearly with propagation distance z (at large 
z). In Figure 22(b) we plot the expansion of the beam width for three sample realiza¬ 
tions of the random waveguide for the weakly disordered case, as explained above. The 
tendency to localize can be seen in each of these samples, but one can immediately see 
the random nature of the process, where different random realizations of the coupled 
waveguide system result in different rates of initial expansion and final beam localiza¬ 
tion width. Lastly, the expansion of the beam width for three sample realizations of 
the random waveguide for the strongly disordered case is shown in Figure 22(c). Note 
the different vertical scale in these three subfigures. It is clear that despite the random 
variation the disorder slows and eventually halts the expansion, and a stronger disorder 
results in a smaller eventual localization width on average. 

It is interesting to study the beam expansion in the Fourier k-space (momentum 
space) as well. A well-known characteristic of the ballistic expansion of the beam is 
that the beam width in the k-space remains invariant under propagation, as shown in 
Figure 23(a). The k-space expansion of the cases plotted in Figures 22(b), (c) are shown 
in Figures 23(b), (c). Interestingly, the beam widths in the k-space conform well to our 
intuition that narrower beams in x-space are wider in k-space', of course, this simplistic 
intuition is not always true, especially in the presence of spatial chirp. 

As was noted above, the statistical nature of Anderson localization means that the 
beam expansion rate and the eventual localization width varies depending on the partic- 
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beam width beam width beam width beam width 



FIGURE 22. This figure shows the expansion of the x-space beam width in the 
(a) disorder-free lattice, (b) weakly disordered lattice, and (c) strongly disordered 
lattice. Three sample realization of the random waveguide are shown for each 
disorder level. The disorder-free lattice shows ballistic expansion, while disorder- 
induced localization is apparent in disordered samples. 



propagation distance (xlOO) 


FIGURE 23. Same as Figure 22, except the beam width is calculated in the k- 
space. The k-space beam width for diffractive propagation in the disorder-free 
periodic lattice of (a) remains unchanged. 



FIGURE 24. Same as Figure 22, except the x-space beam width is averaged over 
100 independent statistical realizations of the disordered waveguides. The error 
bars signify the one standard deviation for the beam width over the 100 samples. 



propagation distance (xlOO) 


FIGURE 25. Same as Figure 23, except the k-space beam width is averaged over 
100 independent statistical realizations of the disordered waveguides. The error 
bars signify the one standard deviation for the beam width over the 100 samples. 
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ular “random” realization of the disordered waveguide. Therefore, no two propagations 
will be identical. However, some general conclusions can still be drawn, in a statistical 
sense, about the average beam width, as well as the variation around the average. In 
Figure 24 we show the average beam width over 100 independent statistical realiza¬ 
tions of the disordered waveguide, plotted as a function of the longitudinal propaga¬ 
tion coordinate z. Figure 24(a) is the disorder-free ballistic propagation and is plotted 
again for comparison with Figures 24(b) and (c), which correspond to the weakly and 
strongly disordered cases, respectively. The vertical scales in these subfigures are all 
chosen identically, for easier comparison. It is clear that, on average, the localized 
beam width is narrower for the case of a stronger disorder, as expected. The error bars 
in Figures 24(b) and (c) signify one standard deviation for the beam width around the 
mean over the 100 independent random samples. Strong disorder results, not only in 
a smaller beam width, but also in a smaller variation around the average beam width. 
Therefore, a stronger disorder is equivalent to a better predictability: when the disorder 
is strong, the beam width is smaller and is almost the same in all random realizations 
of the disordered waveguide. 

Finally, in Figure 25 we plot the corresponding average beam width in k-space, 
where the error bars are now defined as the one standard deviations of the beam width 
in k-space. The interested reader is urged to consult Refs. [37, 53-55] for a more de¬ 
tailed account of the beam expansion in ID and 2D disordered lattices e.g., it is argued 
that in a 2D disordered waveguide the expansion starts as ballistic then goes through a 
diffusive phase, and eventually the localization takes over. However, in ID the diffusive 
transport regime is absent and the expansion turns from ballistic directly to localized. 


Highlights: 

• The light propagating through a coupled array of identical optical 
waveguides remains confined to only a few waveguides (transversely 
localized), if the waveguide-to-waveguide couplings are random. 

• Similar localization behavior can be observed if the individual 
waveguides are randomized, regardless of whether the waveguide- 
to-waveguide couplings are random or not. 

• The stronger the disorder, the more localized the beam is. 


6. Transverse Anderson localization of light in two transverse dimensions 

Transverse Anderson localization in two transverse dimensions can be explored in a 
similar fashion to ID transverse localization discussed in the previous section. Con¬ 
sider the 2D array of optical fibers placed in a hexagonal lattice in Figure 26. The 
hexagonal lattice is special because each fiber has six nearest neighbors separated by 
an equal distance A, and the next-to-nearest separations are \/3A. In contrast, a fiber 
in a square lattice has four nearest neighbors and the next-to-nearest separations are 
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y/2A. Therefore, the next-to-nearest couplings are weaker in a hexagonal lattice and 
the simplifying assumption of the nearest neighbor coupling is more reliable. Similar 




FIGURE 26. A coupled array of optical fibers are placed on a hexagonal lattice. 

Actual simulations in this section are performed on a hexagonal lattice of N = 817 
fibers. 

to Eq. 10, the coupled mode equations for the propagation of the optical field through 
this optical fiber lattice can be expressed as 



(19) 


k 6 NN 


where the sum is only on the nearest neighbors, and the coupling matrix is symmetric. 

As we discussed in the previous section, the transverse Anderson localization can be 
explored by studying the distributions of the elements of eigenvectors of the B matrix 
over the element position space. Here, B is the effective propagation constant matrix 
of the coupled fiber system and is defined as By = /3o Sjj + cjk, noting that cjk are non¬ 
zero only when j and k are the nearest neighbor fibers. The fiber coupling constants are 
assumed to be random in the general form of cjk = co + rjk, where co is a fixed value 
of rjk and belongs to a random distribution. Here, we take j8o = 6 and co = 0.1 on a 
hexagonal lattice of N = 817 fibers. 

In Figure 27 we consider the non-random deterministic case with rjk = 0. The mag¬ 
nitude of the four eigenvectors of the B matrix with the largest eigenvalues are shown 
as density plots on the element-position domain of the hexagonal lattice. For the reader 
familiar with coupled mode theory, the calculated eigenvectors are the supermodes of 
the entire waveguide structure [56], The left most figure shows an azimuthally sym¬ 
metric distribution of non-zero elements that decrease monotonically from the center 
of the lattice toward the boundary. Other eigenvectors have different distributions and 
symmetry properties, but they all nearly fill the entire element-position domain, i.e., 
the non-zero elements of the eigenvectors spread over nearly the entire lattice. 

Fet’s consider what happens when the coupling constants are randomized according 
to rjk G unif[—0.1,0.1] in Figure 28. Similar to the above, the magnitude of four eigen¬ 
vectors of the B matrix are shown as density plots on the element-position domain of 
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FIGURE 27. Each subfigure shows the density plot of the elements of an eigenvec¬ 
tor of the effective propagation constant matrix B on the element-position domain 
of the hexagonal lattice. The density plot shows only the absolute value of the el¬ 
ements of each eigenvector. In the language of coupled mode theory, each density 
plot signifies the intensity distribution of a supermode of the entire coupled fiber 
lattice. This Figure relates to the non-random deterministic situation, where all 
fibers are identical and all coupling strengths to the nearest neighbors are equal. 
Note that the eigenvectors spread over nearly the entire lattice. 



r 











FIGURE 28. This figure is similar to Figure 27, except the nearest neighbor cou¬ 
pling strengths are randomized (off-diagonal disorder). Note that each eigenvector 
is localized in a certain region on the lattice, as expected from the transverse An¬ 
derson localization. 





FIGURE 29. This figure is similar to Figure 27, except the propagation constants 
of the individual fibers are randomized (diagonal disorder). Note that each eigen¬ 
vector is localized in a certain region on the lattice similar to Figure 28, as ex¬ 
pected from the transverse Anderson localization. 
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the hexagonal lattice. All figures show that the non-zero elements of the eigenvectors 
are localized in certain regions on the lattice, as expected from the transverse Anderson 
localization. Of course, the localization region for each eigenvector is different. 

As mentioned before, random variations of the diagonal elements of the B matrix 
can also induce localization of the eigenvectors in the element-position domain. This is 
shown in Figure 29, where the off-diagonal elements are assumed to be constant values 
of co = 0.1 and the diagonal elements (propagation constants of the individual optical 
fibers) are randomized according to B,-,- = j8o + S /i,. where 8 /3, G unif[—0.2,0.2]. There¬ 
fore, the effect of the diagonal disorder is similar to that of the off-diagonal disorders, 
at least qualitatively. Of course, transverse Anderson localization is similarly observed 
if both diagonal and off-diagonal disorder are implemented simultaneously. Transverse 
Anderson localization in a 2D lattice of optical fibers with random variations of the 
parameters was first proposed by S. S. Abdullaev and F. Kh. Abdullaev of the Heat 
Physics Department of Uzbekistan Academy of Sciences in 1980 [29]. 

Earlier in this tutorial review, we briefly highlighted the subtle differences between 
diagonal and off-diagonal disorders in a ID disordered lattice. The weaker localization 
in the middle of the band for off-diagonally disordered lattices persists in 2D as well. 
For example, it has been argued in Ref. [57] that in the middle of the band for a 2D 
disordered square lattice, the wave amplitude falls off with distance as R ' A where R 
is the distance from the center of the mode and A is an exponent that depends on the 
amount of disorder. The power-law localization is also attributed to the geometry of 
the square lattice (and the logarithmic singularity of the density of states in the middle 
of the band). Therefore, for the triangular lattice that has been studied in this section 
(which does not have a logarithmic singularity of the density of states in the middle of 
the band), all modes are exponentially localized. 


Highlights: 

• Transverse Anderson localization is observed for a 2D disordered 
coupled waveguide array, just as in ID. 


7. Transverse Anderson localization of light: the RLV mechanism 

The 2D transverse Anderson localization of light that was proposed by Abdullaev and 
Abdullaev and discussed in the previous section strongly resembles the original pro¬ 
posal by P. W. Anderson. We will refer to this as the AA mechanism for the rest of this 
article, not to be confused with the Aubry-Andre model [58]. In either case, the starting 
point is a lattice-periodic potential where each site on the lattice is characterized by a 
“bound-state” energy (j3o for fibers) and off-diagonal couplings {cjk for fibers). In the 
absence of randomness, the solutions to the wave equations are Bloch-periodic solu¬ 
tions that extend over the entire lattice. If sufficient randomness is introduced in the 
energy and/or couplings, the solutions of the wave equation no longer extend over the 
entire lattice; rather, they localize to certain regions on the lattice. 
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An alternative method to obtain the 2D transverse Anderson localization of light was 
independently proposed by de Raedt, Lagendijk, and de Vries in 1989 [30]. Recall that 
for the AA mechanism, the disorder was introduced on top of an existing ordered lat¬ 
tice. However, in the suggested method by de Raedt, Lagendijk, and de Vries (RLV 
mechanism), the randomness is not superimposed on top of an existing ordered lat¬ 
tice; rather, the underlying potential is completely random [30]. De Raedt et al. [30] 
proposed that this scheme can be realized in a quasi-2D optical system in a dielectric 
with a transversely random and longitudinally uniform refractive index profile. This is 
sketched in Figure 30a, where a long dielectric “waveguide” has a random refractive 
index profile that is invariant along the fiber. De Raedt et al. suggested a binary random 
system, where each pixel is randomly chosen to have a refractive index of n\ or m with 
equal probabilities, where n\ and 112 are marked by red and blue color, respectively in 
Figure 30a. 

Using extensive numerical simulations, De Raedt et al. showed that for a properly 
designed random fiber e.g. with ri\ = 1.0, n 2 = 1.5; and d ~ A. where d is the width 
of each pixel and A is the wavelength of light, they obtain quasi-2D localization in the 
transverse plane of the waveguide. An optical field that is launched in the longitudinal 
direction initially expands until it reaches a terminal localization radius, after which the 
freely propagating beam fluctuates around a relatively stable radius. In Figure 30b, a 
slice of the intensity distribution in the x-z plane is shown for a representative localized 
propagation through a random binary waveguide, where x is a transverse coordinate 
and z is the longitudinal coordinate. 



FIGURE 30. (a) Sketch of the transversely random and longitudinally invariant 
dielectric medium for the observation of the transverse Anderson localization, (b) 
Cross section of a Gaussian beam that is coupled to the disordered waveguide. The 
intensity distribution shows that the beam goes through an initial expansion and 
eventually localizes to a stable width, as expected from the transverse Anderson 
localization. 
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Highlights: 


• One way to observe transverse Anderson localization is to use a 2D 
disordered coupled waveguide array. We refer to this as the AA mech¬ 
anism because it was originally proposed by Abdullaev and Abdul¬ 
laev. 

• Alternatively, transverse Anderson localization can be observed in a 
optical fiber-like medium with a fully random refractive index pro¬ 
file. We refer to this as the RLV mechanism because it was originally 
proposed by de Raedt, Lagendijk, and de Vries. The RLV mechanism 
is distinct from the AA mechanism by the absence the underlying lat¬ 
tice. 


8. Experimental observation of localization: the AA mechanism 

The first experimental observation of the transverse Anderson localization of light was 
carried out in Segev’s group in 2007 [59]. The experiment was more in line with the 
AA mechanism discussed earlier in section 6, where the disorder was superimposed on 
top of an existing triangular lattice of waveguides. The ordered lattice and the super¬ 
imposed disorder were formed, by means of the optical induction technique [60], in a 
photo-refractive crystal (SBN:60) using a laser beam at 514 nm wavelength. 

In order to write the underlying ordered lattice in the photo-refractive crystal, 
Schwartz et al. [59] used the interference pattern generated by three symmetrically 
arranged lasers and obtained a fully periodic hexagonal interference pattern with pe¬ 
riodicity of 11.2 pm inside the 10 mm-long crystal. The disorder is generated from a 
speckled beam that is invariant in the longitudinal direction, as required for the ob¬ 
servation of transverse Anderson localization. The speckled beam was formed from a 
Bessel beam, created by passing a Gaussian laser beam through an axicon, and then 
passed through a 4-f imaging system, where a diffuser was placed at the joint Fourier 
plane between the two lenses. The 4-f imaging system with the diffuser transformed 
the Bessel beam into a broad and speckled beam that was coherently combined with 
the original lattice-forming beam to create the desired refractive index fluctuations on 
top of the ordered lattice in the photo-refractive crystal. 

In order to investigate the localization, Schwartz et al. used another probe beam 
at 514 nm wavelength. The beam width was 10.5 pm full-width at half-maximum 
(FWHM) and was always launched into the crystal at the same location. In the ab¬ 
sence of disorder, they observed diffraction patterns after 10 mm of propagation. They 
observed clear transverse localization of the probe beam when the disorder level was 
increased to more than 30%. In order to obtain an appropriate ensemble to calculate the 
mean localization radius (localization length), the diffuser was rotated by a step-motor, 
such that the laser beam passed through a different location on the diffuser in each step, 
where a new intensity measurement was taken. 

For the large disorder level of 45%, Schwartz et al. fitted the properly averaged inten¬ 
sity distribution of the localized beam, over 100 independent experiment, to an expo- 
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nential of the form exp(— 2\r\/E,), and obtained 15 =64 pm for the localization length of 
the Anderson localized beam. The experiment carried out by Schwartz et al. was quite 
interesting, as they were able to vary the disorder level by controlling the intensity 
level of the disorder-inducing lasers, and also obtained an ensemble of independent 
measurements for the statistical analysis of the localization phenomenon by rotating 
the diffuser. However, the variations of the refractive index of random sites in Ref. [59] 
were on the order of 10 4 . As we will discuss later in greater detail, the radius of the 
localized beam (localization length) depends on the magnitude of the refractive index 
fluctuations. In order to obtain a more localized beam for device applications e.g., com¬ 
parable to the beam radius in a conventional optical fiber, it is necessary to increase the 
magnitude of the index fluctuations [61]. Moreover, when the magnitude of the in¬ 
dex fluctuation is larger, the sample-to-sample variation in the beam radius becomes 
smaller; therefore, each element of the ensemble closely resembles the average. This 
self-averaging behavior alleviates the need for averaging over a large ensemble, and 
for a sufficiently large magnitude of the index fluctuation the average is almost iden¬ 
tical to each element. We will discuss a recent realization of the transverse Anderson 
localization in the presence of a large index contrast in the next section. 


Highlights: 

• The first observation of transverse Anderson localization was re¬ 
ported by Schwartz et al in an AA-mechanism setup. 

• The underlying lattice and the overlaying disorder were optically in¬ 
duced in a photo-refractive crystal. A separate beam was used to 
probe the localization behavior. 

• The setup allowed them to vary the amount of disorder on demand. 


9. Experimental observation of localization: the RLV mechanism 

In section 7 we discussed the RLV mechanism for the transverse Anderson localization 
of light. The quasi-2D system proposed by De Raedt et al. [30] strongly resembles an 
optical fiber. In 2012, Karbasi et al. [17] designed and fabricated an optical fiber that 
functioned based on the RLV mechanism, i.e., transverse Anderson localized in the 
presence of an entirely random refractive index profile. In order to obtain the random 
and pixelated refractive index profile of the RLV mechanism, Karbasi et al. [17] used 
40,000 pieces of a low index polymer fiber polymethyl methacrylate (PMMA) with 
refractive index of 1.49 and 40,000 pieces of a high index polymer fiber polystyrene 
(PS) with refractive index of 1.59. Each fiber was 8 inches long with an approximate 
diameter of 200 pm. The fibers were randomly mixed, assembled to a square preform as 
shown in Figure 31, fused together, and redrawn to a fiber with a nearly square profile 
and approximate side width of 250 pm [62], 

We recall that one of the requirements for transverse Anderson localization is the 
longitudinal invariance. Some of the randomly mixed optical fibers may have crossed 


30 




FIGURE 31. Random mixture of the PS and PMMA fiber strands. 


over each other during the assembly and redraw process; however, the large draw ratio 
of the fiber guarantees that the refractive index profile remains relatively unchanged 
along the fiber. Moreover, small perturbations are likely not going to disturb the trans¬ 
verse localization noticeably. 

Figure 32 shows the scanning electron microscope (SEM) image of a polished poly¬ 
mer Anderson localized optical fiber (pALOF). Figure 32a is the image of the cross 
section of the optical fiber with an approximate side width of 250 pm, where the high 
and low index regions are not distinguishable in this figure. Figure 32b is a zoomed-in 
SEM image of a 24 pm wide region on the tip of pAFOF exposed to 70% ethanol sol¬ 
vent to to dissolve the PMMA, so that the PMMA region can be differentiated by the 
darker color. The sizes of the random features (pixels) in Figure 32b are around 0.9 pm. 



FIGURE 32. (a) cross section of pALOF with a nearly square profile and an ap¬ 
proximate side width of 250 pm; (b) zoomed-in SEM image of a 24 pm wide 
region on the tip of pALOF exposed to a solvent to differentiate between PMMA 
and PS polymer components, where feature sizes are around 0.9 pm and darker 
regions are PMMA; and (c) experimental measurement of the near-held intensity 
profile of the localized beam after 60 cm of propagation through pALOF. The total 
side width of subfigure c is 250 pm, so it can be directly compared with subfig¬ 
ure a. Adapted with permission, copyright 2012, Optical Society of America [17]. 

In order to investigate the guidance and localization properties of the pALOF, light 
from a single mode optical fiber was directly launched (nearly butt-coupled) into the 
pALOF. The output near-field image was collected by a 40X objective and was pro¬ 
jected onto a CCD camera. Figure 32c shows an experimental measurement of the 
near-field intensity profile of the localized beam after 60 cm of propagation through 
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pALOF, where the wavelength of light is 633 nm. The localization was observed to 
be strong: when the input beam was scanned across the input facet, the output beam 
clearly followed the transverse position of the incoming beam [63], 

We encourage the interested reader to consult Ref. [64] for another method to obtain 
transverse Anderson localization based on the RLV mechanism. The optically induced 
randomized potential is created by a computer controlled spatial light modulator (SLM) 
in a photorefractive crystal. 


Highlights: 


• Experimental observation of transverse Anderson localization for the 
RLV mechanism was reported by Karbasi et al in a disordered optical 
fiber medium. 

• The disordered fiber was drawn from a random preform. The preform 
consisted of 80,000 strands of optical fibers with different refractive 
indexes that were randomly mixed and fused together. 

• When the input beam was scanned across the input facet, the output 
beam clearly followed the transverse position of the incoming beam 


10. Detailed analysis of the RLV localization scheme 

An attractive feature of the pALOF is the large index contrast of 0.1 that is helpful 
in reducing the localization radius (localization length) of the beam. The measured 
value for the localization radius was reported as q avg ~31 pat 633 nm wavelength. 
The measured radius was only slightly larger than the calculated value; however, the 
measured standard deviation of a> ~14 pm was considerably larger than the variations 
calculated numerically and reported in Figure 3 of Ref. [17]. The measured standard 
deviation for the beam radius was taken from 100 separate measurements of the beam 
profile radius with 20 different fiber samples and 5 different locations across the fiber. 
This discrepancy between the theoretical and experimental variation of the beam radius 
was mainly attributed to the quality of the polishing of the pALOF tip as the surface 
quality of a polished polymer optical fiber tip is generally lower than what is routinely 
achievable for glass optical fibers. 

Numerical modeling of wave propagation in a disordered fiber 

The simulation of light propagation in a disordered fiber is carried out by numerically 
solving the wave propagation equation Eq. 20 using the finite difference beam propa¬ 
gation method (FD-BPM) [59, 63, 65]. 



( 20 ) 


Eq. 20 is the paraxial approximation to the Flelmholtz equation, where A(r) is 
the slowly-varying envelope of the primarily transverse electric field E(r,t) = 
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Re [A(r)exp(/Ho&oz — i(Ot)} centered around frequency co and ko = 2n/X. n(x,y ) is the 
(random) refractive index of the optical fiber, which is a function of the transverse coor¬ 
dinates, and no is average refractive index of the fiber. The forward propagation scheme 
is implemented using the fourth order Runge-Kutta method [66] . 

The stability condition for the fourth order Runge-Kutta method limits the size of the 
steps in the longitudinal direction as dz < anokodx 2 , where dx = dy is assumed to be 
the size of the transverse grid in the finite-difference numerical scheme, and a = 1 f\p2 
in a uniform medium. In the simulations in Ref. [63], a = 0.02 is chosen in Ref. [63] to 
ensure stability and no power dissipation for reliable long distance propagation. While 
transparent boundary condition [67] is implemented, the size of the simulation domain 
is taken to be large enough to ensure that the total power in the simulation region 
remains unchanged along the fiber for Anderson localized beams. 

A typical reliable simulation of transverse Anderson localization e.g. for pALOF, 
requires a transverse area in the range of ~ 10 5 A 2 . The random refractive index pixel 
size is on the order of a wavelength, and each pixel must be resolved by roughly 10 
points in each direction in the finite difference scheme; therefore, ~ 10 7 points are 
required in the transverse domain with dx ~ A/10. This results in the fourth order 
Runge-Kutta stability criterion of the form dz < aA/10, which is dz ~ A/500 for a = 
0.02. For A «500 nm, dz ~ 1 nm, and for a typical propagation distance of 1 cm required 
for reliable localization, ~ 10 7 steps in the longitudinal direction are required, which is 
computationally intensive. 

The effective beam radius (localization length) is calculated by the variance method 
[30] as 

Uz) = sJ{A( r)|(R-R)2|A(r)>, (21) 

where the angle brackets denote integration over transverse x — y coordinates. R = 
(x,y) is the transverse position vector and R is the vector pointing to the center of the 
beam, defined as the mean intensity position by R = (A(r)|R|A(r)). The optical field 
is assumed to be normalized according to (A(r)|A(r)) = 1/2. 

We note that an estimate of the localization length is sometimes given using a quan¬ 
tity called the Inverse Participation Ratio (IPR) [59]. IPR is defined as 


J I 2 (x,y,z) dxdy 
(f I(x,y;z ) dxdy) 2 ' 


( 22 ) 


where I is the optical intensity defined in the x — y transverse plane, at the longitudinal 
coordinate z. IPR has units of inverse area and the average effective width is defined 
as m e ff = \P) 1 , where (■■•) represents the statistical averaging over the ensemble. 
We note that IPR as defined by Eq. 22 is originally rooted in the definition of the mode 
effective area in nonlinear optics (see e.g. [68]). Equation 21 used by De Raedt et 
al. [30] is likely a better representation of the localization phenomena and more true to 
the random-walk nature of the scattered wave. Therefore, we prefer to use the second 
moment method of Eq. 21, rather than the Inverse Participation Ratio to calculate the 
localization length. 
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Small localized beam radii with low variations are desired 

The experimental realization of the RLV scheme by Karbasi et cii. [17] showed that the 
transversely disordered optical fibers can be embraced as a completely new class of 
optical fibers that guide light, not in a conventional core-cladding setting, but by means 
of the Anderson localization, where any location across the transverse profile of the 
fiber can be used to guide light. As we will see later in image transport applications of 
this disordered optical fiber, a small beam localization radius is likely its main desired 
attribute. A reduced variance is also desired, because as a device it is not helpful to have 
variations in the beam radius that depend on the transverse location of the beam in a 
random and unpredictable way. In other words, although the underlying waveguiding 
mechanism is based on disorder and randomness, the main measurable device attribute, 
the localized beam radius, should be predictable. These considerations lead one to ask 
what can be done, in the design of these disordered fiber in order to minimize the 
localization radius and also reduce the sample-to-sample variations in the localization 
radius. 

Design parameters upon which the beam radius depends 

The localized beam radius can depend on: the input characteristics of the optical beam 
such as the wavelength and the in-coupling beam diameter; the dimensionful fiber char¬ 
acteristics such as the transverse dimensions of the fiber, the width of each pixel; and 
the dimensionless fiber characteristics such as the refractive index of the components, 
especially the index difference, and the geometrical ratio of the components used in the 
fiber characterized by the fill-fraction. 

Given the intuition we established based on the modal picture, we learned that if the 
initial excitation width is smaller than the radius of a typical mode of the disordered 
waveguide (localization radius or localization length) , the propagating light eventu¬ 
ally localizes to a radius that is comparable to the localization length; and if the initial 
excitation is wider than the typical width of an individual mode, the propagating light 
eventually localizes to a slightly wider beam than the original size. Therefore, the de¬ 
sign program to optimize the disordered waveguide for best localization must be inde¬ 
pendent of the input beam profile and must focus on the transverse dimensions of the 
localized beams [69-71], 

Another important property of the disordered fiber is the transverse size of the fiber. 
If the fiber cross section is not large enough, the boundaries of the fiber strongly affect 
the localization radius of the modes. The impact of the boundary on the localization 
radius of the modes that reside near the boundary is inevitable, but in an ideal case the 
fiber cross section must be large enough so that the modes in the interior region are 
shielded from the boundary to prevent excessive scattering loss; therefore, the physical 
characteristics of the majority or nearly all of the interior modes are solely determined 
by the transverse Anderson localization mechanism. 

In the ideal case, when the disordered fiber cross section is much larger than the light 
wavelength, feature size (width of each index pixel), and the localization radius of the 
typical modes, the total transverse size of the fiber becomes irrelevant. In that case, the 
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only two dimensionful fiber design parameters are the light wavelength and the width 
of each pixel. Considering the fact that Maxwell’s equations are scale invariant, if the 
wavelength and feature size (fluctuating index pixel size) are scaled by some factor, the 
localized beam radius of each mode is scaled by the same factor. 

The feature size of the pALOF studied in Refs. [17,63] are approximately 0.9 pm. 
Experimental measurements and numerical simulations have confirmed that the local¬ 
ization radius is much smaller at 405 nm wavelength than at 633 nm wavelength, as 
shown in Figure 33. It is speculated that the optimum feature size for strong localiza¬ 
tion using the RLV scheme is around 2A. This statement is merely an educated guess 
and the optimum value of the feature size relative to the wavelength has yet to be found. 
Intuitively, if the feature size is too small compared with the wavelength, the optical 
field would hop over the random sites and merely average the refractive index, with 
very weak scattering. On the other hand, if the feature size is too large compared with 
the wavelength, the mean free path for wave scattering will be large, resulting in a large 
localization radius. Therefore, there must be an optimum value of the feature size for a 
given wavelength. Finding the optimum value is computationally very challenging. 



propagation distance(mm) 

FIGURE 33. The localized beam radius is smaller at 405 nm wavelength com¬ 
pared with 633 nm wavelength, for the pALOF design with an approximate 0.9 pm 
feature size. Adapted with permission, copyright 2012, Optical Society of Amer¬ 
ica [63], 

Another design parameter that affects the localized beam radius in the binary disor¬ 
dered RFV scheme is the portion of each random refractive index component, charac¬ 
terized by the fill-fraction p. In Ref. [63], p is defined as the portion of the low-index 
material in the higher index host medium. Using an intuitively plausible argument that 
the maximum transverse scattering is obtained when there is an equal amount of low- 
index and high-index material in the disordered matrix, and also a few instances of 
simulations relevant to the experimental regime of interest, Karbasi et al. [63] have 
argued that p = 50% should be regarded as the ideal design target. A sample simula¬ 
tion is shown in Figure 34, where the evolution of the effective beam radius versus the 
propagation distance is shown for different values of the fill-fraction of p = 40% and 
p = 50%; and the latter provides a lower effective beam radius and localization length. 
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propagation distance(mm) 


FIGURE 34. The evolution of the effective beam radius is plotted versus the prop¬ 
agation distance is shown for different values of the fill-fraction of p = 40%, and 
p = 50%, where the latter provides a lower effective beam radius and localiza¬ 
tion length. Adapted with permission, copyright 2012, Optical Society of Amer¬ 
ica [63], 


The last, but definitely not the least important design parameter to reduce the local¬ 
ization radius is the index difference between the random components of the disordered 
fiber. It is generally believed that a stronger localization (smaller average mode radius) 
is obtained if the refractive index difference between the random constituents is in¬ 
creased. This is confirmed in the ID model of the transverse Anderson localization of 
Ref. [46]. However, it is also shown that the dependence of the average mode radius on 
the index difference asymptotically saturates beyond a certain value of the index dif¬ 
ference; therefore, the payback in the reduction of the localization length may be quite 
small beyond a certain threshold index difference. While these issues are well studied 
in the ID model of Ref. [46], the 2D analysis remains to be done. For the 2D geome¬ 
try, Karbasi et al. [63] compared the localization radius in a glass-air disordered fiber 
with pALOF. The refractive index of sites in the glass-air disordered fiber are randomly 
picked as ti 2 = 1.5 and n\ = 1.0 for the glass host and random air-hole sites, respec¬ 
tively. The index difference of 0.5 in the glass-air disordered fiber compared with 0.1 in 
pALOF results in a considerable reduction in the localization radius. In Figure 35, the 
calculated beam radius versus propagation distance for different values of fill-fraction 
arc plotted, and thelocalized beam radius is clearly smaller than that provided by the 
index difference of 0.1 for pALOF. 

Finally, we have repeatedly observed that whenever the scattering strength is in¬ 
creased to reduce the localization length (localized beam radius), whether by selecting 
the right wavelength, or by bringing the fill-factor close to 50%, or by increasing the 
index difference, the statistical variation of the beam radius decreases. Therefore, we 
have the luxury of solving two problems with one solution (a more animal friendly 
expression than killing two birds with one stone!). 
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propagation distance(mm) 


FIGURE 35. Effective beam radius vs propagation distance for different values of 
fill-fraction, p, in glass disordered optical fibers with random air holes. The index 
difference of 0.5 between the random sites results in a very small localization ra¬ 
dius. Adapted with permission, copyright 2012, Optical Society of America [63], 


Spatial beam multiplexing 

The possibility of spatial beam multiplexing in pALOF was studied in Ref. [18] in 
2013. It was shown, both numerically and experimentally, that a pALOF can be used 
to simultaneously transport multiple beams [72]. In Figure 36(a) beam multiplexing in 
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FIGURE 36. Multiple-beam propagation in a 5 cm-long pALOF (a) simulation for 
five beams; (b) experiment for two beams; and (c) experiment for two beams with 
different wavelengths. All beams are at 405 nm wavelength, except the bottom- 
middle beam in subfigure (c), which is at 633nm wavelength. Adapted with per¬ 
mission, copyright 2013, Optical Society of America [18]. 

pALOF is studied numerically at 405 nm wavelength, where the intensity profile of a 
multiplexed beam is plotted after 5 cm of propagation along the fiber; the four exterior 
beams are launched at a distance of 70 pm from the central beam. The output beams 
do not show any appreciable drift and remain in the same transverse location across the 
fiber. 

Experimental verification of beam multiplexing are shown in Figure 36(b), where 
both beams are 405 nm wavelength, and in Figure 36(c), where the upper beam is at 
405 nm wavelength and the lower beam is at 633 nm wavelength. The localization is 
clearly stronger at 405 nm due to the choice of the design parameters of pAFOF, as 
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discussed before. 


Macro-bending loss 

In section 10 the possibility of beam multiplexing in an disordered fiber was discussed. 
Given the unconventional nature of the beam confinement in a disordered optical fiber, 
one must worry about the possibility that the spatially multiplexed beams drift across 
the fiber when the fiber is subjected to substantial macro-bending. This issue was stud¬ 
ied in detail in Ref. [18]. Using numerical simulations, it was shown that transverse 
Anderson localization is very robust and can withstand, at least in theory, a degree of 
macro-bending that is beyond what is acceptable for conventional fibers. To explore the 
macro-bending experimentally, a 10 cm section of a 15 cm-long pALOF was wrapped 
16 times around a mandrel with an approximate radius of 1mm, where no appreciable 
loss of walk-off effect was observed in the localized beam. 

11. Image transport through the disordered fiber 



FIGURE 37. (left) Elements of a group on 1951 U.S. Air Force test target (1951- 
AFTT). Transported images of different numbers through a disordered optical 
fiber: (middle) subfigures (a)-(d) are related to the group 3 on the test target and 
(right) subfigures (a)-(d) are related to the group 5 on the test target. Adapted with 
permission, copyright 2014, Nature Communications [19], 

Motivated by the successful demonstration of beam multiplexing, Karbasi et al. used 
pALOF for endoscopic fiber-optic imaging. To their pleasant surprise, the image trans¬ 
port quality was comparable to or better than some of the best commercially available 
multicore imaging fibers, with less pixelation and higher contrast [19]. Figure 37 shows 
some of the transported images in the form of numbers from a section of the 1951 U.S. 
Air Force resolution test chart through pALOF. The test-target, in the form of a stencil 
in which numbers and lines were carved, was butt-coupled to the hand-polished input 
facet of pALOF and was illuminated by white light. The near-field output was projected 
onto a CCD camera with a 40 x microscope objective. 

The minimum resolution of the images is determined by the width of the point spread 
function of the disordered optical fiber imaging, which was calculated to be smaller 
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than 10 jam at 405 nm wavelength [63]. In practice, the imaging resolution in pALOF 
is limited by the quality of the cleave and polishing of the fiber. The fiber surface quality 
is partially responsible for the distortions in the transported images in Figure 37. The 
high quality image transport in the proof-of-concept experiment has been achieved 
without any optimization in the design. 

Highly multicore optical fibers, similar to those introduced in section 6, have been 
used for direct transportation of images in various configurations [73-76]. The highly 
multicore fiber, often referred to as the “coherent fiber bundle,” is commonly used in 
medical and industrial endoscopy [73,77,78]. 

For an Anderson localized fiber, a higher amount of disorder and a larger level of 
fluctuation in the refractive index provides stronger beam localization, resulting in an 
improved image resolution. A similar statement can be said for multicore imaging fiber 
as well: the coherent core-to-core coupling is detrimental and blurs the image; there¬ 
fore, core-to-core coupling must be suppressed by varying the size of the cores so 
that neighboring cores cannot couple resonantly. Also, in a multicore imaging fiber, 
the cores must be as close as possible to create a less pixelated image, while being 
very different in size so that they cannot couple efficiently even when they are close. 
One can view the disordered Anderson localized fiber as taking the two limits of the 
high-packing of the cells and the large variation of the sizes to an extreme, so that the 
individuality of each core is completely lost: all neighboring sites are strongly coupled, 
but the extreme randomness prevents light leakage and blurring. 



FIGURE 38. Transported images through the disordered fiber and the commercial 
image fibers. Images related to group 5 of the 1951-AFTT test chart in (a) pALOF, 

(b) FIGH-10-350S image fiber and (c) FIGH-10-500N image fiber (experimental 
measurements). The scale bar in (b) is 30 /Im-long and the same scale bar can 
be used for (a) and (c). Each fiber is approximately 5 cm long. Adapted with 
permission, copyright 2014, Nature Communications [19]. 

The imaging performance of the “unoptimized” pALOF compares with some of the 
best commercially available multicore imaging optical fibers, as is further confirmed 
in Figure 38. The transported images over 5 cm of the number “6” from group 5 of 
the 1951-AFTT test chart are compared between pALOF in Figure 38(a), Fujikura 
FIGH-10-350S in Figure 38(b), and Fujikura FIGH-10-500N in Figure 38(c). The im¬ 
age quality of the transported image through the pALOF is clearly better than FIGH- 
10-350S and is comparable with FIGH-10-500N. The feature sizes in Figure 38 are on 
the order of 10-20 pm. The Rayleigh range for this level of resolution is approximately 
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0.8-3 mm, which is substantially shorter than the typical propagation length in these 
imaging fibers. Therefore, the imaging results are non-trivial and cannot be obtained 
using bulk propagation or conventional multimode fibers. 


Highlights: 

• In an Anderson localized disordered optical fiber, a smaller localized 
beam radius is obtained via increasing the differences between the 
refractive indexes of the random dielectric constituents. 

• The feature size should be ~ 2A, and a fill-fraction of 50% is pre¬ 
ferred. 

• A stronger localization is generally accompanied by a smaller vari¬ 
ation in the localized beam radius, resulting in a more uniform and 
predictable beam radius. 

• Beam multiplexing and high-quality image transport are featured as 
device-level applications of the transverse Anderson localization in a 
disordered optical fiber. 


12. Transverse Anderson localization in a disordered silica optical fiber 

The first observation of Anderson localization in a silica fiber was reported in Ref. [79]. 
The main motivation for using a glass-air structure has been the larger index contrast 
that results in a smaller beam diameter (better image transport resolution), as well as a 
lower sample-to-sample variation in the value of the beam diameter (better image uni¬ 
formity) [79]. The reported glass-air disordered fiber was drawn at Clemson University. 
The preform was made from “satin quartz” (Heraeus Quartz), which is a porous artisan 
glass. By drawing the preform, the airholes (bubbles) in the glass are stretched to form 
the hollow air-rods required for transverse Anderson localization. The large draw ratio 
sufficiently preserves the longitudinal invariance, without significant disturbance over 
typical lengths used in the experiments. 



FIGURE 39. (a) SEM image of the glass optical fiber with random airholes re¬ 
ported in Ref. [79]; and (b) zoomed-in SEM image of the same fiber. 

The cross-sectional SEM image of the disordered glass-air optical fiber is shown in 
Figure 39a, and a zoomed-in SEM image is shown in Figure 39b. The SEM images 
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provide a good estimate of the refractive index profile of the fiber; the light gray back¬ 
ground matrix is glass and the black dots represent the random airholes. The diameter 
of the disordered glass fiber is measured to be 250 pm. The diameters of the airholes 
vary between 0.2 pm and 5.5 pm. Unfortunately, the airhole fill-fraction was shown to 
be as low as 2% in the central regions, which was far below the ideal value of 50%, so 
strong localization was not observed near the center. However, Anderson localization 
was observed near the boundary, where the airhole fill-fraction is 8%. 

More recently, transverse Anderson localization has been reported by scientists from 
Corning Incorporated in random airline fibers (RALF) [80]. The preform for RALF is 
fabricated using the outside vapor deposition (OVD) process. The consolidation pro¬ 
cess is done in the presence of 100% nitrogen gass, during which nitrogen is trapped 
in the blank to form glass with randomly distributed air bubbles. When the preform 
is drawn, the random air bubbles stretch to form random airlines. RALFs with 150, 
250 and 350 pm diameters were fabricated, where the averaged air line diameters were 
177, 247 and 387 nm in these fibers, respectively. The maximum airhole fill-fractions 
measured in RALFs is reported to be 1.473%, which is considerably lower than that 
of Ref. [79]. Ref. [80] reports the observation of transverse Anderson localization of 
light, despite the low airhole fill-fraction, and attributes its observation to the substan¬ 
tially smaller diameter of the airholes compared with that reported in Ref. [79]. 

We also encourage the interested reader to consult a pioneering work by Pertsch et 
al [81] on light propagation in a disordered 2D array of mutually coupled optical fibers. 
They observed both localized and delocalized modes and also analyzed the impact of 
nonlinearity on these modes. 


Highlights: 

• Glass-air random optical fibers are highly desired because the large 
refractive index difference between glass and air can result in nar¬ 
rowly localized beams, and also small variations around the mean 
beam radius. 

• All attempts at glass-air random fibers have so far resulted in unde¬ 
sirably small air fill-fractions. 


13. Hyper-transport in longitudinally varying disordered waveguides 

As mentioned in section 10, the paraxial propagation of an optical beam in a longitu¬ 
dinally invariant medium is described by the paraxial approximation to the Helmholtz 
equation, Eq. 20. This equation is formally equivalent to the Schrodinger equation in 
quantum mechanics, where the time variable t is replaced with the longitudinal coor¬ 
dinate z and the potential is replaced by the term proportional to the refractive index 
profile n(x,y) 2 — The longitudinal invariance of the refractive index profile is taken 
into account by explicitly showing that n is only a function of the transverse coordi¬ 
nates x and y. The requirement for the longitudinal invariance of the refractive index 


41 




profile is equivalent to the potential being constant in time in the Schrodinger equation. 

Presence of some form of time dependence in the random potential can potentially 
change the localization picture. This issue has been studied extensively over the years 
e.g., in Refs. [82-86] for temporally random or correlated potentials. An interesting 
finding has been the possibility of hyper-transport in certain temporally varying ran¬ 
dom potentials, where the wavefunction expands even more rapidly than the ballistic 
expansion observed in free space. We recall our earlier discussions on ballistic, dif¬ 
fusive, and localized propagation of light in disordered optical waveguides. For the 
propagation of the optical beam, the beam width generally grows with the propagation 
distance according to w(z) °< z p , where p = 1 for ballistic, p = 1/2 for diffusive, and 
p = 0 for localized propagation. An example of the ballistic expansion is the familiar 
formula for the width of a Gaussian beam as a function of the propagation distance [39] 

w(z) = „,,y , i + (|) , a = (23) 

where w(z) z for z zo- 

For a temporally random potential, Ref. [83] showed analytically that the width of 
the wavefunction, as measured by the second-moment method, expands with time ac¬ 
cording to f 3 / 2 . Hyper-expansion was later shown for correlated temporally fluctuating 
potentials as well in Refs. [84-86]. 

Hyper-transport for an optical beam was studied experimentally and theoretically 
in Ref. [87] and Ref. [88], respectively, where the temporal variations were replaced 
by rapid longitudinal fluctuations in the refractive index profile of the waveguide. The 
experiment reported in Ref. [87] was performed in a similar setting to that of the first 
observation of the AA mechanism earlier reported in section 8 and Ref. [59]. The dis¬ 
ordered lattice was formed by means of the optical induction technique [60] in a photo- 
refractive crystal. However, unlike Ref. [59], where special care was taken to ensure 
the longitudinal invariance of the interference speckles and the resulting transversely 
random index profile, the refractive index profile created in the nonlinear crystal in 
Ref. [87] was made to vary with z. The experiment proved that in the presence of a suf¬ 
ficiently large longitudinal variation of the index fluctuation, not only did the transverse 
localization cease to exist, but also the beam expanded at a rate faster than ballistic as 
it propagated through the random-index waveguide. 


Highlights: 

• In the presence of rapid longitudinal fluctuations in the refractive in¬ 
dex profile of a disordered waveguide, it is possible for the optical 
beam to expand even more rapidly than the ballistic expansion ob¬ 
served in free space. 
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14. Anderson localization and nonlinearity 

The inteiplay between nonlinearity and disorder is of great interest to the studies of 
transverse Anderson localization. The main question that is often asked is whether 
the presence of nonlinearity perseveres, enhances, undermines, or destroys Anderson 
localization. This issue has been explored over the years, where some of the studies 
have benefited from the existing literature on systems with similar dynamical equations, 
such as the Bose-Einstein condensate in the presence of disorder (see e.g. Ref. [89] and 
the references therein). Here, we briefly highlight some of the results that are more 
relevant to the concept of transverse Anderson localization. For more details, we refer 
the interested reader to an excellent review on this subject by Fishman, Krivolapov, and 
Soffer, in Ref. [90]. 

Earlier, in section 8, we reviewed the numerical and experimental work of Schwartz 
et al. [59] that resulted in the observation of transverse Anderson localization of light 
for the AA mechanism. The authors also investigated the transverse Anderson local¬ 
ization of light in the presence of Kerr nonlinearity, both numerically and experimen¬ 
tally. The defining equation for the nonlinear propagation of light is the nonlinear 
Schrodinger equation (NESE) and is identical to Eq. 24, with the addition of a Kerr 
nonlineaity term [68]: 

dA 1 

i ——|- ——— [Vj-A + £5 («“ — 715) A] + ko»2|A|“A = 0 , ( 24 ) 

dz 2«oko 

where «2 is the nonlinear index, which is positive for self-focusing and negative for 
self-defocusing nonlinearity. 

As a case study, the authors considered a disordered lattice where the maximum 
contribution of the nonlinear term to the index change “max (|/?2 x |A| 2 )” was assumed 
to be a maximum of 15% of the index contrast of the underlying periodic waveguide. 
They also varied the disorder level from 0% to 30%, where the disorder level was de¬ 
fined as the magnitude of random index fluctuations relative to the index contrast of the 
underlying periodic waveguide. They observed that over this range, the self-defocusing 
nonlinearity ( 722 ) results in a moderate (nearly negligible) widening of the average beam 
profile. However, the self-focusing nonlinearity (n 2 > 0) resulted in a substantial reduc¬ 
tion of the average localized beam diameter. The enhancement of localization due to 
the self-focusing nonlinearity was particularly noticeable when the disorder level was 
less than 15%. 

The experiments were carried out at 15% disorder level. However, the maximum 
nonlinear contribution was taken to be equal or higher (up to a factor of 3), compared 
with the index contrast of the underlying periodic waveguide. This was achieved by 
making the probe intensity equal to or higher than the interference maxima of the 
lattice-writing beams (for a perfect lattice). The statistical analysis of the localized 
beam radius clearly confirmed the expected reduction in the average beam radius due 
to the self-focusing nonlinearity. 

Similar results were reported by Lahini et al. [53] using disordered one-dimensional 
waveguide lattices. Their - experiment consisted of a one-dimensional lattice of coupled 
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optical waveguides patterned on an AlGaAs substrate. Light was injected into one or 
a few waveguides at the input, and light intensity distribution was measured at the 
output. We recall our earlier discussions in section 5: we identified the highly localized 
eigenmodes near the top edge of the propagation constant band. The amplitude of these 
flat-phased modes, as they are referred to in Ref. [53], are in-phase at all sites; perhaps 
only with a few flips as seen in Figures 20(a) and (b). Another set of highly localized 
modes were shown to exist near the top edge of the propagation constant band. The 
amplitude of these staggered modes, as they are referred to in Ref. [53], have phase 
flips between adjacent sites, as shown in Figures 20(e) and (f). 

In the weak nonlinear regime, Lahini et al observed that nonlinearity enhances local¬ 
ization in flat-phased modes and induces delocalization in the staggered modes. This 
behavior is explained as follows: the presence of the weak nonlinearity perturbatively 
shifts (increases) the value of the propagation constant of each localized mode. For the 
flat-phased modes, the nonlinearity shifts the modes outside the original linear spec¬ 
trum. However, for the staggered, which belong to the bottom edge of the propagation 
constant band, a perturbative increase in the value of the propagation constant shifts 
it further inside the original linear spectrum. Therefore, the propagation constant of a 
staggered mode can cross and resonantly couple with other modes of the lattice, result¬ 
ing in delocalization [91-94]. 

In 2008, Pikovsky and Shepelyansky presented a somewhat different account of the 
interaction between disorder and nonlinearity [89]. Here we rephrase their main find¬ 
ings in a language more consistent with our notation so far. In a disordered coupled 
waveguide lattice, they demonstrated that above a certain critical strength of nonlinear¬ 
ity the Anderson localization is destroyed and turns into a subdiffusive spreading. They 
focused on the discrete Anderson nonlinear Schrodinger equation, which is essentially 
the same as Eq. 10 with the addition of a third-order diagonal Kerr nonlinear term, 
expressed as 



Pj ) A j(z) + c o A j+ 1 (z) + A j- t(z) 


+ y\Aj(z)\ 2 Aj(z)=0, 


j = h -,N. 

(25) 


They also assumed that the waveguide coupling coefficients co is deterministic and is 
identical for ah waveguides. The disorder is introduced through the diagonal propa¬ 
gation constant terms, where they are assumed to be randomly distributed according 
to 


Pj F unif[j3 0 - y, P° + y ] • (26) 

For the linear case of y = 0, the modes are exponentially localized due to the disorder. In 
the nonlinear case where y / 0, Pikovsky and Shepelyansky demonstrated that above 
a certain critical strength of nonlinearity, the Anderson localization is destroyed and 
the field spreads in a subdiffusive form indefinitely across the optical lattice. They also 
showed that the mode width calculated using the second moment method (see Eq. 21) 
grows with the propagation distance as z a , where 0.3 < a < 0.4. 
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Pikovsky and Shepelyansky presented a theoretical argument in support of the subd- 
iffusive spreading of the beam with a = 0.4 [95]; however, their work is mainly based 
on numerical integration of Eq. 25 and monitoring the results, up to z = 10 8 /co- For 
the numerical simulation they used the boundary condition Aj(z = 0) = 8jj m , where 
j m represents the middle waveguide, and the integration is performed by the operator 
splitting method. They also assumed that Y!J=i |Ay(z)| 2 = 1, without any loss of gener¬ 
ality. In a sample set of simulations they chose the nonlinearity strength to be y = co 
for two cases: case 1 with SB = 2co; and case 2 with SB = 4co- The second moment was 
evaluated according to 



As expected, the initial expansion was ballistic for either case, but after some distance 
Zo, the expansion became subdiffusive. They tit the subdiffusive expansion to a(z) = 
(Joz a over the range z,o < z. < 10 8 /co- In case 1, for different instances of randomness, 
they obtained 0.32 < a < 0.39; and for case 2 they reported 0.28 < a < 0.41. Upon 
averaging over 8 independent realizations, they reported a tit of the form 57.5 x z 0,344 
for case 1 and 8.7 x z 0306 for case 2 over the subdiffusive range. They also reported 
a critical value of nonlinearity y c ~ O.lco above which this subdiffusive behavior is 
observed. 

Intuitively speaking, one may think that in a nonlinear disordered coupled waveg¬ 
uide system, the dynamics of the beam is initially influenced by nonlinearity; and as 
the beam spreads, the effect of nonlinearity becomes weaker and the disorder dynamics 
takes over. Therefore, one should always expect Anderson localization after sufficiently 
long propagation. This is clearly in contrast with the findings of Pikovsky and Shep¬ 
elyansky reported above. Fishman et al. [90] present a thorough survey of the many 
subtleties involved rearding the interaction of nonlinearity and disorder. The conclu¬ 
sion is that the situation can best be described as inconclusive at this point. For ex¬ 
ample, when using the numerical simulations, they caution that Eq. 25 is chaotic with 
an exponential sensitivity to numerical errors. For long-distance propagation, it is not 
clear that reducing the z step size can control the cumulative numerical error, given that 
the limit of zero z step may be singular. Details are beyond our intended scope and can 
best be found in Ref. [90]. 

Other forms on nonlinearity besides Ken - can also interact with the disorder-induced 
localization. For example, it was recently shown that a beam of light propagating in 
a pAFOF (introduced in section 9) exhibits self-focusing properties due to a thermal 
nonlinearity [96]. The larger light absorption strength in PMMA than PS results in a 
inhomogeneous temperature distribution. The higher temperature in PMMA translates 
into a decrease of its refractive index. The result is an increased refractive index mis¬ 
match and stronger localization. The results are quite counter-intuitive, because the 
polymer materials used in the experiment have defocusing intrinsic nonlinear coeffi¬ 
cients (112 < 0). In Ref, [96], Feonetti et al. demonstrated that transversally localized 
modes shrink when the pump intensity is increased despite the fact that 112 < 0 for the 
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polymers. In a subsequent publication [97], the authors provided further evidence of 
this behavior by analyzing the direct relation between the optical intensity and the lo¬ 
calization length, and also demonstrated the disorder-induced focusing by a monochro¬ 
matic continuous wave (CW) laser. 

The interested reader is also urged to study other aspects of the interaction between 
nonlinearity and disorder not covered in this tutorial review e.g., on soliton propagation 
in random media [98,99]. 


Highlights: 

• The main question of interest is whether the presence of nonlinearity 
perseveres, enhances, undermines, or destroys Anderson localization. 

• The answer to the above question can best be described as inconclu¬ 
sive at this point. 


15. Coherence, classical and quantum light, and Anderson co-localization 

In the previous sections Anderson localization was explored for a temporally and spa¬ 
tially coherent light, coupled to a disordered optical lattice. In this section some of the 
main issues related to the propagation of partially incoherent light are explored. A dis¬ 
cussion on the propagation of single photons and correlated photons is also presented. 

Anderson localization of waves with imperfect coherence was reported in 2011 by 
Capeta, et al [100]. The main issue that was answered in their paper was the extent to 
which transverse Anderson localization is affected by the partial coherence of the in- 
coupling beam in a disordered linear lattice. They observed that if all the eigenmodes 
of the disordered waveguide are exponentially localized, any partially incoherent beam 
exhibits localization with exponentially decaying tails, after sufficiently long propaga¬ 
tion distances. The reported observation conforms with intuition, because an incoherent 
wave can be thought of as a superposition of coherent modes with stochastically vary¬ 
ing coefficients. Because each coherent mode is expected to undergo localization, the 
entire beam should localize as well. However, localization is delayed by incoherence 
compared with the case of a coherent in-coupling beam: the more incoherent the wave 
is, the longer it diffusively spreads while propagating in the medium. 

For an optical beam propagating in a ID disordered optical lattice described by the 
propagation Eqs. 10 or 20, the state of the coherence of the beam is determined by the 
mutual coherence function r (1 )(xi,X 2 ,z) = {A* (x 2 ,z)A(x\ 1 z)) [39], where (•■■) is the 
ensemble average, and A is the stochastic held. For the disordered waveguide array of 
Figure 17, x\ and xi are the waveguide array indexes. For example, the mutual coher¬ 
ence function of the partially coherent extension of the input Gaussian beam of Eq. 18 
can be written as 



exp ,_UzAf^zAi ]exp| . 


U~k ) 2 


j,k= ,N, (28) 
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where jo is the index of the waveguide at the center of the input Gaussian beam; and 
and "A) are the spatial and the coherence widths of the beam, respectively. We note that 
the optical intensity over the coupled waveguide array is given by the diagonal element 
of the mutual coherence function F'-'j (j is the waveguide index) [39]. 

Using Eq. 10 and following a similar procedure that resulted in Eq. 17c, we conclude 
that the propagation of the mutual coherence function is described by 

r$(z) =it bf^W^k^M-KPf-Mz], (29) 

k '=t /=i 

where the propagation coefficients are determined at z = 0 according to 

bj',k' =tt r 5*(z = 0)vf } . (30) 

k'—l /—1 

Let’s consider the disordered coupled waveguide array in Figure 18, for /3o = 6, 
Co = 0.01, N = 201, rj 6 unif[—0.1,0.1], and 0 < z < 15000, according to the notation 
used in the discussion following Eq. 10. The input light is assumed to be a partially 
coherent beam defined by Eq. 28, with Wq = \Jl and jo = 101. The average beam 
width as a function of the propagation distance is plotted in Figure 40 for the case of 
highly coherent .Aq = 100 in red, semi-coherent ,Aq = 5 in blue, and near-incoherent 
.A() = 2 in cyan. The beam widths are averaged over 100 independent simulations. The 
results agree with the observations of Capeta, et al [100], where localization happens 
for all three cases, but is strongest for the more coherent input beam. 



FIGURE 40. The average beam width as a function of the propagation distance in 
a disordered coupled waveguide array is plotted in for the case of highly coherent 
,A) = 100 in red, semi-coherent -A) = 5 in blue, and near-incoherent /A) = 2 in 
cyan. Wq = \/2 and jo = 101 have been used and the beam widths are averaged 
over 100 independent simulations. The disordered coupled waveguide array is 
defined by /3 q = 6, cq = 0.01, N = 201, rj £ unif[—0.1,0.1], 


The propagation of non-classical light in disordered waveguides has also been a sub¬ 
ject of interest over the past few years. For example, the propagation of a Fock number 
state, coherent input state, as well as a squeezed state in a disordered waveguide array 
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was studied by Thompson et al [14]. The quantum mechanical analogue of the disor¬ 
dered tight-binding model used for the propagation of light in a disordered waveguide 
array in Eq. 10 can be heuristically constructed by simply elevating the amplitudes of 
the optical held Aj to photon annihilation operators dj in waveguide j, subject to the 
following commutation relations 


[dj,d k \ = 0, 


= 8 kj , j,k=l,---,N. 


(31) 


The quantum mechanical analogue of Eq. 10 is the “linear” Heisenberg equation for 
the evolution of dj(z) and dj(z) and the creation and annihilation operators at the output 
ports can be obtained from those at the input ports using the Green’s function of the 
disordered medium as 


aj(z) = Z G jk(z)d k ( 0). (32) 

k= 1 

The Green’s function, in the language of the disordered couple waveguide array and 
Eq. 17c, can be expressed as 


G jk (z) = £ Vf Vf } exp[-ift*]. (33) 

k'=1 

Using Eq. 32, the input first-order quantum coherence function can be mapped to 
that at the output according to 


r$(z) = I £ (G'jkWMz)) d W rg,( 0 ), (34) 

k— \ k'—\ 

where r[^,(z) = (aj(z)d£,(z)). The averaging on the Green’s functions are earned over 
multiple realizations of disordered waveguides hence the subscript “dw” in (• • • )d w . 

Equation 34 is exactly what one would write for the propagation of the first-order 
classical coherence function as well. The quantum effects enter through the first order 
coherence function, which is defined by the density matrix p as 

r S'(°)= Tr {p«1(0)^'(0)}. (35) 

For example, if the input state to each waveguide is a Glauber coherent state, and the 
quantum state of the input light is defined by | a \, 0 , 2 , ■ ■ • , a#), we will have 

E$(0 ) = a* k a k ,. (36) 

Alternatively, if the input state consists of Fock number states defined by 
\ni,ri 2 , • • • ,«;v), we will have 


r$(0 ) = n k s kk >. 


(37) 
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The average output intensity in port j, Ij{z) is the diagonal element of the output 
first-order coherence function. For the case of the coherent state above, we will have 


N N 


Ij(z) = £ £ (G* k (z)G jk ,(z ))dw a* k a k . 


(38) 


k=lk'=1 


This formula is identical to the intensity we would expect to obtain at port j when 
classical optical field of amplitude a k are coupled at the input port k. For the Fock 
number states, the intensity is given by 


N 


h(z) = Yj ( G *jk(z)G jk (z)) dw n k . 


(39) 


k= 1 


The special case where the photons are coupled into a single input port, say middle port 
M, results from Eqs. 38 and 39 can be simplified as 


Ij(z) = {\GjM(z)\ 2 )dw |«m| 2 , or = (\GjM(z)\ 2 )dw nut- (40) 


Equation 40 is notable-it shows that regardless of the statistics of the input pho¬ 
tons, if it is injected to a single port we will observe the same localization behavior. In 
fact, localization is solely dictated by the coherent averaging of the Green’s function 
of the disordered waveguide (|G/,v/(z)| 2 )dw This is not the case when light is coupled 
into multiple ports, as can be observed in the difference between Eq. 38 and Eq. 39: the 
coherent interference of the input coherent states in Eq. 38 certainly affects the localiza¬ 
tion, while Eq. 39 merely consists of multiple copies of Eq. 40 assembled incoherently 
together because of the total ambiguity in the phase information in the Fock number 
states. The bottom line is that the localization is dictated by both the input photon statis¬ 
tics as well as the disordered waveguides, and their interaction can be complicated or 
simple depending on the situation. 

In Figure 41 we plot the \Gj k (z)\ 2 at z = 10,000 for the disordered coupled waveg¬ 
uide array in Figure 18, for /3o = 6, co = 0.01, N = 201, rj E unif[ — 0.1,0.1]. We also 
assume that k = 101; therefore, the input light is assumed to be coupled only to the 
middle waveguide as discussed in Eq. 40. \Gj k (z)\ 2 in Figure 41 is plotted as a func¬ 
tion of j, which labels the output waveguide number. The dashed blue line represents a 
single simulation, while the solid red line shows the result of averaging ( \GjM(z)\ 2 )dw 
over 1000 simulations. The exponential localization is clearly observed in Figure 41. 

Higher order quantum coherence functions can also be studied in disordered waveg¬ 
uides. In particular, the second-order quantum coherence function gives us information 
about two photon correlations and Hanbury Brown-Twiss (HBT) effect [101]. The 
second-order quantum coherence function at the output of the disordered waveguide 
can be written as 



(41) 


N N 


= I I (G* k (z)G* fk ,(z)G fV (z)Gj l (z)) dw (4(0)a+(0)a>(0)d z (0)) qs . 


k=iw =t 
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FIGURE 41. \Gjk[z)\ 2 for k = 101 is plotted as a function of the output waveg¬ 
uide number for a disordered coupled waveguide array. The dashed blue line rep¬ 
resents a single simulations, while the solid red line shows the result of averaging 
(|G/m(z)| 2 ) dw over 1000 simulations. 


We are interested in two particular scenarios. The first scenario is when the two 
photons are coupled to the same waveguide, say middle port M, and the input state 
is given by (d^) 2 |vac). In this case, the second-order quantum coherence function is 
given by 

I?>) = (\G jM (z)\ 2 \G fM (z)\ 2 ) dw . (42) 

The second scenario is when the two photons are coupled to two different waveguides, 
say ports M and M', and the input state is given by a^ M ,a' M \ vac). In this case, the second- 
order quantum coherence function is given by 

rj.‘l(z) = (\GjM(z)Gj>M'(z)+ GjM'(z)GfM(z)\~}dw (43) 

In either scenario, we have indistinguishable photons that co-propagate in the lattice, 
and are subject to the “averaged” correlations of the propagators of the disordered 
lattice, resulting in rather exotic quantum statistics behavior, as shown for example 
in Ref. [13]. In the following, we will adopt the same disordered waveguide array as 
that of Figure 41 and illustrate the behavior of the second-order quantum coherence 
function. 

( 2 ) 

In Figure 42 we plot T\ j, (z) of Eq. 42 for the case when the two photons are coupled 
to the middle waveguide M = 51 (we use N = 101 waveguides for easier simulation), 
in the absence of any disorder (rj = 0), where z = 2400 and z = 10,000 is assumed in 
Figures 42(a) and (b), respectively. The plots indicate the ballistic expansion of the pho¬ 
tons, with larger probability to find the photons in the edge waveguides (Figure 42(a)), 
and also reflections from the boundary for longer propagation (Figures 42(b)), resulting 
in strong interference effects. 

In Figure 43 we present plots of the same scenario of Figure 43(a), except in the 
presence of disorder, averaged over 1000 independent simulations. We have used o-e 
unif[—0.002,0.002] for Figure 43(a) and rj € unif[—0.004,0.004] for Figure 43(b). We 
have also presented both density and 3D plots in each case for easier comparison. Note 
that we have taken z = 2400 in order to prevent reflections from the outer waveguides; 
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FIGURE 42. rH(z) of Eq. 42 is plotted for M = 51 (N = 101 waveguides in the 
lattice) as a function of the output waveguide numbers j and / for a disorder-free 
coupled waveguide array for (a) z = 2400, and (b) z = 10,000. 


such reflections complicate the analysis and are not observed as long as there is a 
sufficiently large number of coupled waveguides in the array for a given propagation 
distance. 

In the presence of low disorder as in Figure 43(a), it is possible for both photons to 
remain localized (center of the plot), one remains localized and one propagates freely 
(edge-centers of the plot), or both freely propagate (corners of the plot), in decreasing 
order of probability. Once the disorder is increased as in Figure 43(b), it is almost 
only possible to have localized photons with small non-zero chance of having one 
localized and one freely propagating photon. Once the disorder is further increased e.g., 
rj £ unif[—0.01,0.01] (not plotted here), the only reasonable non-zero probability is to 
have both photons highly localized. We recall that the disorder-free case in Figure 42(a) 
favors both freely propagating to one localized and one freely propagating to both 
localized, respectively, and the presence of disorder reverses this behavior. 

In Figure 44 we plot Y ~ (z) of Eq. 43 for the case where the photons are coupled to 
different waveguides: M = 50 and M' = 51 in Figure 44(a) and M = 50 and M' = 52 in 
Figure 44(b). Either case is in the absence of any disorder (rj = 0) for z = 2400. The 
plots indicate the ballistic expansion of the photons, with larger probability to find the 
photons in the edge waveguides (Figure 42(a)), and also reflections from the boundary 
for longer propagation (Figures 42(b)), resulting in strong interference effects. 

When the two photons are coupled to neighboring waveguides as in Figure 44(a), 
the most likely scenario is that the two photons ballistically propagate to the same edge 
of the lattice, and it is impossible to find the two photons in the opposite edges of 
the lattice; however, there is a small yet non-vanishing probability to find one photon 
remaining in the center waveguide and the other in the edge lattice. When the two in- 
coupled photons are separated by one waveguide as in Figure 44(b), the most likely 
scenario is that the two photons ballistically propagate to either edge of the lattice, and 
it is impossible to find one photon remaining in the center waveguide and the other in 
the edge. 
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FIGURE 43. r (2) ,(z = 2400) of Eq. 42 is plotted for M = 51 (N = 101 waveg¬ 
uides in the lattice) as a function of the output waveguide numbers j and / for 
a disordered coupled waveguide array for (a) r ; - £ unif[—0.002,0.002] and (b) 
rj £ unif[—0.004,0.004]. Both density and 3D plots are presented in each case for 
easier comparison. 


In Figure 45 we present plots of the same scenario of Figure 44(a), except in the 
presence of disorder, averaged over 1000 independent simulations. We have used 
rj £ unif [—0.002,0.002] for Figure 45(a) and rj £ unif [—0.004,0.004] for Figure 45(b). 
We have also presented both density and 3D plots in each case for easier comparison. 
We note that increasing the level of disorder compared with the disorder-free case in 
Figure 44(a) decreases the probability of ballistic co-propagation of photons to the 
same edge, while increasing the probability of the localization of both photons in the 
center waveguide. Therefore, the disorder-free case of Figure 44(a) favors both freely 
propagating to the same edge to one localized and one freely propagating to both lo¬ 
calized, respectively, and the presence of disorder reverses this behavior. 

We encourage the interested reader to consult Ref. [13] for a more detailed account 
of the behavior of upon which much of the preceding discussion is based. 

The above discussion only covers the case of off-diagonal disorder in a disordered 
coupled array waveguide. The comparison with the case of the diagonal disorder is 
performed in Ref. [15]. It is shown that the evolution and localization of the photon 
density (optical intensity) is similar in the two cases of diagonal and off-diagonal dis¬ 
order, as expected. However, the intensity correlation (photon density-density corre- 
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FIGURE 44. r (2) , (z) of Eq. 43 is plotted for (a) M = 50, M' = 51 (N = 101 waveg¬ 
uides in the lattice), and (b) M = 50, M' = 52, as a function of the output waveg¬ 
uide numbers j and / for a disorder-free /V = 101 coupled waveguide array and 
z = 2400. 


( 2 ,) 

lation or P. j, (z)j carries a distinct signature of the type of disorder. The propagation 
of an entangled-photon pairs in a disordered waveguide array has also been explored 
in Ref. [16]. It has been shown that while neither photon is localized, the two-photon 
separation in coincidence space is: this behavior is called Anderson colocalization. 
The increase in entanglement is accompanied by a gradual evolution from Anderson- 
localization to Anderson colocalization. 


Highlights: 

• If all the eigenmodes of a disordered waveguide are exponentially lo¬ 
calized, any partially incoherent beam exhibits localization with ex¬ 
ponentially decaying tails, after a sufficiently long propagation dis¬ 
tances. 

• The more incoherent the in-coupling beam is, the longer it diffusively 
spreads while propagating in the medium. Therefore, localization is 
delayed by incoherence compared with the case of a coherent in¬ 
coupling beam. 

• The observed localization behavior in a disordered waveguide is af¬ 
fected by the quantum statistics of the in-coupling photons. The ef¬ 
fect can be mainly observed in the second-order quantum coherence 
function, which gives us information about two photon correlations. 


16. Conclusions 

A tutorial review of the transverse Anderson localization of light in disordered waveg¬ 
uides is presented. In addition to the intriguing fundamental science behind Anderson 
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FIGURE 45. r , / 2, / (z = 2400) of Eq. 43 is plotted for M = 51 and M' = 52 (N = 101 
waveguides in the lattice) as a function of the output waveguide numbers j and f 
for a disordered coupled waveguide array for (a) rj £ unif[—0.002,0.002] and (b) 
rj £ unif[—0.004,0.004]. Both density and 3D plots are presented in each case for 
easier comparison. 


localization, there are potential applications for light interaction with disordered dielec¬ 
tric systems. Examples such as spatially-multiplexed beam delivery and high-quality 
image transport were explored. There are still plenty of open questions in both linear 
and nonlinear dynamics of disordered systems that have yet to be answered. There has 
been a recent surge of interest among the optics community in Anderson localization, 
and it would have been impossible to cover them all even briefly in this survey. This 
brief tutorial review is intended to demystify Anderson localization for the broader au¬ 
dience and provide sufficient background to newcomers in this field to follow up with 
studying other more in-depth reviews and research articles. 
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